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• PREFACE 


The main results of our analysis of the transport of 
cosmic rays in the interplanetary and interstellar media are 
studied and summarized. We show that previous modulation 
theories are incorrect except for high energies (> 1 BeV) 
where solar modulation effects are observed to be very small. 

The main outcome of our effort is a new convection-diffusion 
equation which is valid in the energy region where solar modu- 
lation is observed. 

We formulate the problem of the transport of cosmic rays 
in full generality from first principles (Chapter II ). A 
master equation is deduced which allows for a closed formula- 
tion of the transport problem. We show that two dimensionless 
parameters (a measure of particle energy and a measure of field 
fluctuation strength) characterize the transport regimes. In 
the ultrahigh energy regime where the Fokker-Planck equation 
holds, contact is established with conventional theory. We 
then resum the set of contributions that corresponds to only 
two-point correlations with straight line trajectories. This 
resummed equation allows us to extrapolate from ultrahigh to 
high energies. 

We then deduce. the equations for the omnidirectional inten- 
sity and flux as the first two spherical harmonic components 
of the velocity distribution function (Chapter III). These 
equations are the basis for transport theory in the modulation 
energy range. To discuss these equations, we develop suitable 
asymptotic expansions (Chapters IV, V) which are applied to 
deduce the cosmic ray transport theory. 

We strongly recommend that our transport theory be applied 
to the wide range of modulation problems of current interest 
since it is the only derived theory which is applicable in the 
modulation region. 

A summary of special problems and publications is given 
on page x immediately preceding the main body of this report. 
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SPECIAL PROBLEMS AND PUBLICATIONS 


This report covers work carried out during the three-year 
period of January 12 , 1968 to January 11, 1971* The main pur- 
pose of our program has been to implement Fermi's proposal that 
many important properties of cosmic rays should be understood 
in terms of the interactions of charged particles with the 
turbulent electric and magnetic fields contained in the large 
scale plasma that surrounds the Sun and penetrates the galactic 
disc. Our work on cosmic ray modulation is discussed thoroughly 
in this report. In addition, we have also been concerned with 
a number of special topics; in particular: 


' 1 - Electrostatic effects in stellar winds. 

2 - The heating of hydrogen clouds by cosmic rays and 

the ‘ possibility of giving a noncosmo logical explan- 
ation for the black-body background radiation. 

3 - Statistical mechanics in relativistic form, including 

as a special case relativistic thermodynamics. This 
-subject needs considerable clarification, especially 
when electric fields play a significant role. Even 
the simplest model • (two-body relativistic potentials) 
presents substantial difficulties. 

4 - A theoretical study of x-rays from the solar wind 

in the low-energy band . 

* \ 

5 - A feasibility analysis of a window less x-ray counter 

whose main chamber is protected by a flow of gas 
through a fluid dynamic nozzle . 

We include below, a list of the publications that have been 
made during the past three years. The main body of the text 
‘following the list of publications discusses our work on the 
transport theory of cosmic rays. 
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CHAPTER I 


INTRODUCTION - OUTLINE OP THE PROBLEM 
AND SUMMARY OP RESULTS 



This report contains our development of the only cosmic 
ray transport theory available which is derived rigorously 

s 

from first principles and is suitable for studies of the solar 
modulation problem. The only assumptions made in our develop- 
ment of this theory are: (l) we have neglected strong anisotropy 
in the cosmic ray beam, (2) we have neglected higher tha-h two- 
point correlations in the random magnetic field, and (3) we have 
assumed the two-point correlation tensor to have an isotropic 
form. 


The result of our analysis can be stated in the form of a 
convection-diffusion .transport equation for the cosmic ray omni- 
directional intensity, with the relationships between the trans- 
port coefficients and the mean and random magnetic fields given. 
We summarize this result here. The transport equation is 

|§ + V. (VI) = ^ V.^-VI (1) 

where I is the omnidirectional intensity which is a function of 
position, time, and rigidity. All lengths in Eq. (l) are 
measured in units -of L , a macroscopic length which measures 
the size of the region in which the particles are contained; 
time is measured in units of L/v , where v is particle speed. 
The diffusion tensor has two components for diffusion -perpen- 
dicular and parallel to' the mean magnetic field; they are given 
by 


and 



( 2 ) 


D 


II 



(3) 


The convection velocity, V , can be shown to be an average Alfven 
drift velocity due to curvature and gradients in the mean field. 
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but modified due to the interaction with the random field. It 
is given by 

v = v X (pv) (4) 

A 

where P is a unit vector in the direction of the mean field and 


V = 



( e+a) g) 2+a> ± 


(5) 


The quantities co^ , and co g are related to the correlation 
function through the solution of the following dispersion rela- 
tions 


A 


to,, 


- -C«n) a /“ ax e '“" 


(6) 


to. 


and 


co 


- -(enff 

J o 


co —tO. X 

6 Xa 


K, (e,X)cos co g X - K g (e,\)sin co^X 


& J 


(7) 


= (en) s f 

on 


00 -to X 
dXe 


JjC_l( £ > x ) sin co g X + K g (e^) cos co g X 


g vlQ 

The correlation function is written in the isotropic form 


( 8 ) 


We define 


and 


R ij (?) - s ij A(r) + B ( r ) r i r j 


*,(£,>■) B (tf) f dzA(r)P^(z) 

= (^W)p <( z ) 


(9) 


(10) 


(ii) 
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where is the i'th-order Legendre polynomial and 


■(«.*.») = x[z 2 +(l-z 2 ) (\~ °°| ? X ) ' 


1/2 


( 12 ) 


We now introduce the following four combinations 


t 4ir 
J 1 ' 3 


_ A 0 " 5 A 2 J * J 2 3 [5 A 2 + A o] ^ 3 ) 


rr _ 4 t |[ 

K 1 “ 15 


sr b 4 ' 7 B s + B 


0] 


(14) 


rr _ 4 t£i 
K 2 ~ 15 1 


21 B 4 + 2 B 2 + B 0 


] 


(15) 


The kernels in Eqs. (6) through (8) are given by 


K H (e,X) = 2J X cos eX - 


( 16 ) 


K ± (e,X) = J 1+ J 2 cos eX-X 2 K 2 2 ( - s - 1 g r ^ --^l- - )-2( xy°. .. y . r~± J-i 


(eX)‘ 


(IT) 


and 


K (e,X) = J _ sin eX 

o 


x 2 k. 


Vcos eX-lV^n sin eX^ 

IA — sc — A 1 ' r - e^yJ 


( 18 ) 


These equations give a transport description of the cosmic ray 
particles in random magnetic fields. The complexity of these 
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equations is necessary to describe the low rigidities (less than 
1 BV) in which modulation effects are observed. For very high 

2 

rigidity, our results simplify considerably to^ = -(eq) 

and (D = 0 
g 

There are two independent parameters in this theory. They 
are: e, a measure of the inverse rigidity of the particles, and 

q, a measure of the strength of the fluctuations in the field. 

Our work is based on a truncated master equatiai for the cosmic 
ray distribution function which we derive from the Liouville equa 
tion using the technique of Kaufman.'*' This equation can be 
written schematically in the form 

§| = (-eq) 2 C(e)f (19) 

where D/dt is a convective time derivative taken along the 
trajectory of the particles in the mean field and C(e) is an 
integrodifferential operator which gives the effect of the random 
field on f .By deriving this equation from first principles, 
we have found that the right side of Eq. (19) represents the 
leading term of an expansion of the exact master equation in 
powers of (eq) . Higher-order terms in (eq) bring in the effect 
of higher-order correlations in the random field . In the inter- 
planetary field, e=l at approximately 1 BV . Since e is 
inversely proportional to rigidity, we see that to discuss solar 
modulation effects we must develop a theoiy which is valid for 
e > i . 

Previously developed theories of cosmic ray transport in 
random magnetic fields can be shown to be based on equations of 
the form of Eq. (19) which are, in some cases, identical to our 
truncated master equation and, in other cases, have been shown 
to be equivalent . 2 > 33^*5 However, in all of these theories, when 
modulation energies are discussed, the approximate effects of 
C(e) are determined by expanding C in powers of l/e and 
retaining the leading term. These expansions are in clear 
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violation of the condition « 1 unless i] -+0 . In the 

transport theory which we have presented in the preceding para- 
graphs, no expansion in e has been made. Thus, the theory is 
valid in the range e ~ 1 

In the remainder of this report, we describe the derivation 
of our transport theory. In Chapter II, we formulate the problem 
of the transport of cosmic rays in full generality from first 
principles. We keep explicitly only magnetic fields, though 
electric fields could be introduced in most of this chapter with- 
out great difficulty. The problem is set up in terms of the 
Liouville equation which is then ensemble-averaged over the dis- 
tribution of magnetic fields. A master equation which consti- 
tutes the foundation of our work is deduced which allows for a 
closed formulation of the transport problem. Two special cases 
are then discussed. In the ultrahigh energy regime, contact is 
established with conventional theory. We then resum the set of 
contributions that corresponds to only two-point correlations 
with straight line trajectories. This resummed equation has an 
H-theorem and allows us to extrapolate from ultrahigh to high 
energies . 

In Chapter III, we deduce the equations for the omnidirec- 
tional Intensity and flux as the first two spherical harmonic 
components of the velocity distribution function discussed in 
the previous chapter. These equations are the basis for the 
transport theory in the modulation energy range. With the tech- 
niques 1 of the next two chapters, we have deduced the new con- 
vection-diffusion equation given here. 

In Chapter IV, we develop the asymptotic expansion for 
systems of linear equations. The techniques developed in this 
chapter apply to the master equation of Chapter II, as well as 
to the equations of Chapter III. We show how to extract the 
secularities and how to obtain regular expansions. 
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We find that for short-range two-body correlation functions, 
when e ~ 1 , the flux equation is an integrodifferential equa- 
tion with a long-range kernel. Kernels with short range have 
been treated previously and the available uniformization tech- 
niques are adequate. The special problems presented by long- 
range kernels are discussed and solved in Chapter V. 

In Chapter VI, we introduce electric as well as magnetic 
fields. In particular, we obtain complete solutions for the 
problem of charged particles moving in uniform fields. Two 
Lorentz- in variant, gauge- invariant constants of the motion are 
obtained which turn out to' be adiabatic invariants. We, thus, 
have laid the groundwork for an expansion in small electric 
fields (such as those in Alfven waves). 
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CHAPTER II 

GENERAL FORMULATION OF THE PROBLEM - 
ULTRAHIGH AND HIGH-ENERGY APPROXIMATIONS 



1 

Following Fermi's original paper in which he suggested 
that the cosmic ray flux is of galactic origin but still iso- 
tropic due to the presence of a disordered galactic magnetic 
field, a great deal of work has been done to refine our under- 
standing of the transport of energetic particles in random 
electromagnetic fields. In the earliest of this work, transport 
equations were assumed which usually had one feature in common: 
the connection between the relevant transport coefficients and 

O 

the state of the random field was only loosely modeled. Attempts 
to strengthen the connection between the transport coefficients 
and the random field have been made, quite naturally in studies 
of the kinetic theory which might underlie -the cosmic ray trans- 
port theory. In some cases,^the Boltzmann equation has been 
adopted, but then it has become necessary to devise an effective 
cross section for the particle-random field interaction. Even 
so, this approach has led to transport theories in which all the 
transport coefficients could, at least, be related to the assumed 
cross section. A major advance was made when several investiga- 
tors introduced the Fokker-Planck equation as a possible kinetic 
theory.^" In the transport theories which resulted, the transport 
coefficients could be directly related to measurable properties 
of the random field . 

Recently, attempts have been made to justify the assumption 
of the Fokker-Planck equation as the relevant kinetic theory for 
energetic particles in random electromagnetic fields. These 
attempts have been based on the more fundamental Liouville equa- 
tion or its equivalent. Hall and Sturrock^ using the quasi- 
linear approximation, have constructed a kinetic theory from the 
Liouville equation which they have demonstrated to be equivalent 
to the earlier Fokker-Planck theories. They use a test particle 
analysis in which the effects of the moving charged particles on 
the electromagnetic field are neglected. Kulsrud and Pearce^ 
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have developed a more general, self-consistent theory in which 
the field can be affected by the energetic particles. In both 
of these theories, under suitable conditions, the authors have 
shown that they are able to reproduce the results of Jokipii^ 
which are based on the Fokker-Planck equation. The essential 
conditions which are necessary for these comparisons are: (1) a 
stationary random magnetic field with no electric field and 
(2) low rigidity particles so that the radius of gyration of the 
particles is negligibly small compared to the correlation length 

associated with the random field. 

\ 

The work which we present in this report falls Into the class 
of papers in which the appropriate cosmic ray kinetic theory is 
derived from the Liouville equation. Our work does not embody 
all of the generality of the previous papers in this class; we 
'limit ourselves to a test particle- analysis of particles moving 
in a stationary magnetic field with no electric field. Early in 
our work, we reproduce a special case of the kinetic theory of 
Hall and Sturrock; however, our approach is considerably differ- 
ent. We do not use the quasi-linear approximation at the outset; 
and, as a result, we are able to determine more exactly the con- 
ditions under which this kinetic theory is a valid approximation. 
We find these conditions in contradiction to the assumption of 
low rigidity particles which is made by Jokipii, Hall and 
Sturrock, and Kulsrud and Pearce when they attempt to evaluate 
the Fokker-Planck coefficients or their equivalent. This is not 
to say that the Fokker-Planck or the equivalent kinetic theories 
do not apply; the problem lies in further approximations which 
have been made in the evaluations of the Fokker-Planck coeffi- 
cients. 

7 

In our approach, following the technique of Kaufman,' we 
construct from the Liouville equation the exact master equation 
for the ensemble-averaged, one-particle distribution function. 

We then expand the master equation in powers of the magnitude of 
the random field. Stopping at second order, we obtain the 

9 



kinetic theory of Hall and Sturroclc with no electric fields, 
which is equivalent to the Fokker-Planck theory as well. The 
conditions under which the relevant expansion parameter is small 
can be stated as follows. Let us measure the particle radius of 
gyration in the mean field in units of the two-point correlation 
length in the random field. When the radius of gyration is of 
the order of the correlation length, then r ~ 1 . Xf rj is 

O 

the ratio of the root-mean-square value of the random field to 
the magnitude of the mean field, then the Fokker-Planck theory 
is a valid approximation so long as r g » -q . The previous 
attempts of Jokipii, Hall and Sturrock, and Kulsrud and Pearce 
to evaluate the Fokker-Planck coefficients by obtaining the 
leading-order term in an expansion of the coefficients in r . 

o 

about r =0 is clearly in violation of this condition. 

o 

. In Section 1, immediately following this Introduction, we 
derive the master equation. In Section 2, we expand the master 
equation to second order in eq . Also, in Section 2, we 
expand our "interaction operator," the equivalent of the Fokker- 

Planck coefficients, up to n'th order in e and then resum. 

8 9 

We use the time-scale extension technique J ^ to ensure the 
compatibility and uniformity of our expansions; a simple pertur- 
bation expansion fails for large times. In Section 3^ we note 
that an H-theorem can be calculated from our kinetic theory; an 
arbitrary initial distribution function always relaxes to a 
final, isotropic state. In addition, -we calculate the relaxation- 
time necessary for the final state to be reached . From these 
calculations, we find an additional constraint on our theory. 
Distribution functions which are very anisotropic cannot be 
described by it. 
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1 . DERIVATION OF THE MASTER EQUATION 


We ignore the effects of interparticle collisions and, 
instead, concentrate on the behavior of charged particles in a 
magnetic field characterized by a field ensemble mean part, 

, and a random part, B” . Then the Liouville equation 

AM M4 

for the N-body distribution function reduces to 

-1- KF^ + e(£ + = 0 (l.l) 

where F v '(x,p,t) is the one-body distribution function which 

AM M\ 

depends on position, x , direction on the unit-momentum sphere, 

A VtAA 

p , and time t . The differential operators which appear in 

wt 

Eq. (l.l) have all been written in dimensionless forms in the 
following manner. A length, X , has been introduced which is 
a measure of the two-point correlation length in B' . Time is 
measured in units of the traversal time, X/v , where v is the 
particle speed. Thus 

t = t (v/X) ( 1.2) 

Spatial gradients in f are measured in units of X ; therefore 

K= i-Wr (i - 3) 

The Lorentz operators £ and £' are written as 


£ = _p.fi.iL_ 

^ dp 


{!.>) 




and 


£ = -p • 0, 

jiAA AM, 


MA 


•( 1 . 5 ) 


where Q and D. are the skew- symmetric tensors defined by 
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( 1 . 6 ) 


< 3 ^> 

j ” € ijk^k “ € ijk <B>“ 


-and 



e iO‘k P k " € ij'k 


B k 

V <b"'~b'> 


(l.7> 


The carat denotes the unit vector, p = p/p , and, in addition, 
we use the notation tybp ± = p^-^p^/Sp.j). Both £ and z' do not 
operate on functions of the magnitude of the momentum. Dynamical 
variables in momentum reduce to the two directions on the unit- 


momentum sphere. 


As a result of this dimensional analysis, - two parameters 
emerge . The first of these 


T1 


/<b'.b'>/<b> 


( 1 . 8 ) 


is a measure of the relative strength of the random part of the 
field. The second parameter is 


• P 


(1.9) 


where P is the particle rigidity, e can be thought of as the 
ratio of the correlation length to the gyro-radius of the particle 
in the mean field . 

We look for solutions to Eq.- (l.l) of the form 

^ f + g (1.10) 

where f is the field ensemble average of p(^^ and g repre- 
sents the fluctuations in p( 1 ^ . to determine the relevant 

kinetic theory for f , we further assume that f is • independent, 
of position. In a later paper, we will allow weak gradients to 
determine the transport properties of the particles in the random 
magnetic field. 

Following the technique of Kaufman 1 (1968) we have determined 
from Eq. (l.l) that f obeys uhe following master equation 
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( 1 - 11 ) 


efif = (€T 1 ) 2 (l- e r 1 < ( C'G>)“ ;i <£ 

where 

G ^ ^ + K + e£ + eq£ (1.12) 

Thus, G is the Green's function operator with the total 
Hamiltonian 

» = K + + eq£' ( 1 . 13 ) 

for the generator of the particle motion. Before going on, we 
note that Eq. (l.ll) is considerably more general than implied 
here. Any force field for which a Hamiltonian for the particle 
motion can be constructed could be included in the formalism 
which leads to Eq. (i.n). In particular, in subsequent work, 
we plan to include electric fields with random parts. 


2. EXPANSION IN eq WITH LINEAR TIME SCAIES EXTENSION 


Given the master equation for f from the previous section, 
we notice that it contains two independent parameters e and 
eq . In this section, we construct an asymptotic expansion of 
f to second order in eq . We employ the linear time scale 
extension technique to ensure the uniform validity of the expan- 
sion for large time. 

A power series expansion of the master equation in eq is 
given by 

^ + e£f = (eq) 2 <£ Gq£ >f + O(eq)^ (2.1) 

where Gq is the Green's function integral operator which inverts 

°0 1= lc +K+eC -(2-2) 


13 



The generator of the particle motion is H = K + e£ , the 
Hamiltonian operator in the mean field. Notice that the second- 
order interacticn operator involves only the two- 

point correlation function in the random field. Thus, an 
expansion of f to second order in erj neglects the effects 
of higher-order field correlations. 

A generalization of Eq. (2.1) in which electric fields 
have been included has been derived by Hall and Sturrock.^ 

They use the quasi-linear approach to go directly from the 
Liouville equation to their version of Eq. (2.1), bypassing 
the construction of the master equation. They have also pointed 
out the equivalence of the Fokker-Planck equation and Eq. (2.1) 
to second order in erj . Indeed, they find that their theory 
for stationary magnetic fields and in the limit of vanishing 
rigidity reduces to the low rigidity work of Jokipii^ which 
is based on the Fokker- Planck equation. Kulsrud and Pearce 
find that under the same conditions, which they adopt for most 
of their paper, their Fokker-Planck theory of cosmic ray 
propagation reduces also to the low rigidity results of 
Jokipii . 

We find ourselves in disagreement with the low rigidity 
results of these previous theories. By deriving Eq. (2.1) 
from the exact master equation, we find the specific condi- 

Q 

tion for its validity as an approximation, (erj) « 1 
The low rigidity results of Hall and Sturrock, Jokipii, 
and Kulsrud and Pearce are computed asymptotically in the 
limit of vanishing 1/e , where Eq. (2.1) itself is not 

correct. 
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Han ana Sturrock ana KuXa.ua ana 

specific, comparisons ^f tHeir theor^e^ ln tha 

results of xouipia vhrch . ^ results wWh those of 

limit of vanishing term of our e-expansion gives 

jolcipii's, voe find tha i<a in agreement with 

hls high rigidity results Thus,-. ^ rlgialtys ou r 

jokipii's asymptotical j tlons fco jokipii 1 s when e - 1 

theory contains signi.ica ^ oorrelat ion length in the 

(when the radius of gyra i ° of magnitude) and our theory 

rTai™ with the expansions in V* • 

have attempted a tJ. The^ormal 

ii - ~ r"“ rrv“ rr. =;;»=- 

of f a we employ th~ tim . 1T , hation expansion. In the 

this secular behavior from the per perturbation terms will 

following, the nonuniform behavior of .h. pe-tu 

be clearly exhibited. h 

We introduce an extended distribution function,,! 
is a function of many independent time sea e* ’ O^VS ^ 

£ ? c considered an extensile ^ « ^Litidimensional 

f(t t r t s ,..) reduces to *M 

We further expand f in powers of eh 


2 

f * f Q + «(£]. + ( £T| ) iz + 


(2.3) 


and use the freedom we ^ 

in the time-scale space to ensure 

from Eq. (2.1). is uniformly valid in time. 

we have found that a linear trajecto^ is auxtahle^^ ^ 
removing nonuniform terms from the expansi 
The trajectory is given by 
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(2*4) 


and along this trajectory 


Ot 


c) 


(en) 


dx^' 


lr£ + 


(2*5) 


The extension of Eq. (2=1) is obtained by substituting Eqs. (2*3) 
and (2|.5) into Eq. (2.1); we do not extend the interaction oper- 
ator. By combining the coefficients of each power of (srn), we 
find to second order 


s £o 


eZ to 


= 0 


( 2 . 6 * 


*V- 




^ + + e£ -i 


= 0 


(2.7); 


and 




Sf 0 


+ £fif a = <£'G 0 C'>f_ 0 


From Eq. (2.5) we find 


( 2 * 8 ) 



-e£ x> 




(2,9) 


We have used the following notation: f 0 (T 0 ) gives the dependence 

of f Q on Tq and all higher-order time scales to be introduced 
while “ o(^l ) Si V3S the dependence of £q on and all the 

higher-order time scales but with x Q =- 0 . By substituting 

Eq. (2.9) into Eq. (2.7) and integrating over Tq , we find 


f , 5 So<*o> 

-l( T 0^ 0 ’St; 


( 2 . 10 ) 


To remove this secular behavior on the t q scale from f 1 we set 
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and 


£i(*o> ■ 0 


( 2 . 11 ) 



In this case* Eq. (2.8) reduces to 


Sf 2 

^ + .^2 



4- Gq<£ >f Q 


( 2 . 12 ) 


(2.13) 


We have been unable to evaluate the secular behavior of the 
perturbation term in Eq. (2.I3) directly. The difficulty is rela- 
ted to the helical trajectory generated by h . On the other 
hand , a further expansion of G Q in powers of e is possible. 

We find 


G q [ 1 + £ G qq£] 1 G qo - (-e) n (G 0C) £) r G 00 (2-14) 

n=0 

where the generator of the particle motion in Gqq is simply K; 
the particle moves along a straight line at constant momentum. 
This simplification in the particle trajectory is actually suffi- 
cient to overcome the added complexity introduced by the infinite 
summation in Eq. (2.14). Upon expanding f Q and fg in powers 
of e , we find we are able to evaluate the secular behavior 
of each term in this expansion of the interaction operator. We 
find that a nested time-scale extension in e is necessary to 
ensure ihe uniformity of the expansions of f Q and f 2 

Substituting Eq. (2.14) into Eq. (2.13)5 we find 

00 

^ j 

+ e£ £ 2 = + Zj ( G 00 £ ) nG 00 <C >£0 ( 2 . 15 ): 

0 2 n=0 


IT 



where 


<* C ( G 00 £ ) nG 00 ,C ^ f 0 <asK n (s)£ / (x-ps J p)>f 0 (p,s) (2-16) 


with 


« n (s) ’ / dXK n-iW^W 


(2 .,17) 


for n > 0 and with J£ n (s) - 1 . The quantity £ is defined by 


£(t) S e- K ' c £e K ' t 


(2.18) 


which reduces to 


£(t) = £ - x k. 


( 2 . 1 ^). 


where 


A Ji A 


/e~ = p*£2»V = p*f2— t — 

2 dx/x 

We now expand and . fg in powers of e so that 


( 2 . 20 ) 


-0 ~ ^00 + € -01 + 6 -02 + * ‘ * 


( 2 . 21 ) 


and 


-2 ~ -20 + e -21 + € -22 + * • * 


( 2 . 22 ) 


We further extend these functions and choose the trajectory 


/ „\m n 

T mn " ^ s T 


(2.23) 


so that, in particular 5 
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a 


dT 0 ^ 




+ e 


2 d 


01 


ct 


+ . 


02 


(2.24) 


and 


3 _ 3 . c 5 

3~ - Z7Z; + 6 3? 


+ £ 


2 d 


+ . 


k 2 ~ 20 
Prom Eq. (2*6). we find 


21 


22 


(2.25) 


d -00 


OTT 


= 0 


00 


(2.26) 


and then 


df df 

I ~01 ^00 _ = 0 

* T 00 ' O,r oi “ 00 


(2.27) 


Since does not depend on t qo , 


-01 


= - T 


00 


■b£ M 

-00 

l ° t 01 


+ £f 


00 


( 2 . 28 ) 


Notice that we are restricting our choice of initial conditions to 

-0i( T 00 =0 ) = 0 (2.29) 

for all i > 0 «. To remove the secular growth in f^ , we must 

let 


and then 


£22 + £,f 0 

otq^ .*00 


£01 "° 


‘Again,, from Eq. (2.6) 


^-02 + ^-00 
^00 ^02 


- 0 


(2*30) 


(2.31) 


(2.32) 
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Therefore 


(2-33) 


(2.34) 


(2-35) 


(2.36) 

From Eqs. (2.15), (2.1 6) and (2*17), we find to zeroth order 

in € 

3?^ ■ - 37^ + < £ ( X ’P> J 0 < 3 s^s)£ (x-ps, p )>£ q0 (p) (2.37) 

where, from Eq. ( 2.26 ), ?qq does not depend on the variable of 
integration, s . In Appendix A we show that i"f 

■A a ✓ j * 

r “ X - X = -ps .(2.38) 

and if we choose the isotropic form for the two-point correlation 
function for the random field 10 

R i j (? ) = A(r)6^ + B(r) r ± r^ .(2-39) 

where A- and B are arbitrary even functions of r , then 

<£'(x,p)£'(x,'p)> = A(r)V 2 -(2.40) 


f ^-00 

~02 T oo 


We- let 


^-00 

° and f Q2 = 0 
This process can be repeated to obtain 


f Qi = 0 i = 1,2,3' 


and 


Ir 2 = 0 i = 2.3. 

OT 0i 


where 
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(2.41) 


V 


2 


3 

A. 


A 

ap 


Is the angular part of the Laplacian operator in spherical momen^ 
turn coordinates. Thus, we find 



^00 

We assume that A(s) 



is sharply cut off at s 


(2.42) 


1 so that 


f T °° 

/ dss n A(s) (2.43) 

J 0 

exists and is finite for all n as » . On the other 

hand , if 

p T 00 

/ dsA(s) 

J 0 

***** 

reaches a constant finite value for Tqq > 1 , then integration 

of Eq. (2.42) over t q0 to obtain f 2Q will lead to growth of 
the perturbation term directly proportional to t 00 for large 
t 0 q . The time scale extension technique has been introduced 
by us to remove this secular behavior of the perturbation terms. 
We set 


a £oo 


37 


20 


= J 0 


(2.44) 


and 


, , r'-OO e 00 0 

~20^ T 00^ ~ -20^ t 01^ ” J n ds J d£'A(s')v 2 f 00 (2.45) 

u s 


which leads to a finite £ 2 0^ T 00^ ( an(3 > therefore, a uniform 

expansion of f^ to this order in e and further gives us the 
leading kinetic equation for f Q in e . The nonzero £ 2 0^0l) 
must be reatined to allow a uniform expansion in higher orders of 
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To first order in e , we find from Eq. (2.15) 


d -20 Zf b -21 

37 ^ -20 


^-00 


OT 


21 


<£'(x J p) J dsK l (s)£'(x-ps,p)>f 00 (p) (2.46) 


By direct computation, we find that 

[V?£] = 0 

Therefore, from Eqs. (2.45) and (2.30)* 


(2.4 7 ) 


^ 20 ^ 00 ^ 




01 


+ £ - 20 ^ t 00 ^ = 


^-20 ^ t 01^ „ . 

” + ^-20^01^ (2.48) 


01' 


The first-order kernel is 


Kj>) = s£ - i s 2 k 2 


(2.49) 


To evaluate the effect of k 2 in the interaction term, we write ' 




<£'(x,p)fc 2 £'(x'p)> = - p^<n^.(x)/c 2 «' m (x')>p 

^ p i 


% aB 


(2.50) 


m 


Now /c acting on u (x-ps) can be written as — Z acting on 
, A /\ c ~ s 

ft (x-ps) . Therefore 


<£'(x,p)/c 2 <c'(x;p)> = - K- p j [u<^ j .(x)o' iti (x')>! (2.51) 

p i ' p m 

where the symbol j j J signifies that £ acts only on the quan- 
tities within the bracket and not on everything to the right. The 
expectation value can be evaluated and on carrying out the rest 
of the operations indicated in Eq. (2,51) we find 
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Thus 


<£'(x,p)e 2 £'(£.p)> = 0 


(2.52) 


^ 20^01 ^ M ,± \ . ^-21 




01 


+ ££ 20 (t 0 i) + 37 


00 


^-00 , 
^21 


- <£'(x,p) J 00 dss££'(x-ps,p)>f 00 (p) (2.53) 


and 


-21 ^ t 00^ “ 


-21^01^ 


^£20^ T 01^ , ^ ^ 
^ooV^tt: + ^20^01 v 


01 


- T 


00 37 


b £qo 


- r°°ds f ds (x,p)s'££.' (x-ps^p)yf 

On On ■ _ 


21 ^0 


00 


(2.54) 


The procedure at this point is as follows. The perturbation term 
can be written in the general form 


^1 ( t oo^£oo " t 00^ 7 1 +A 1^00 ( 2 *55) 

where Qj(tq 0 ) a bounded function of t 0 q and 7^ and 
are tTqq independent. They have been differentiated by the 
condition 


[vfa^ = 0 
[V®7 X ] / 0 


(2.56) 


We set 


£21^00) “ £21^01^ Q 1 ^ t 00 ^£ 


■21 1 01' 


1 ^ 00 y -00 


(2.57).- 
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(2.58) 


-2CA T 01 

^oi 


+ £ -20^ T 01' 


7 l-00 


and 

^-00 

37^ - A l£oO ( 2 *59) 

This separation of terms is dictated by the condition that 

,Ji . Q i 

-21 T OCr be bounded and ^ tk® compatibility condition 


d 2 f 

-00 

d-r^dT 


00 ~ 21 


d 2 f 


37 


“°37t: - t^£oo - 0 


21 00 


x (2.60) 


The secular part of the perturbation term which commutes with 
~2 

V determines the kinetic behavior of f Q0 on the Tg^ scale , 
while the secular but noncommuting part determines £20^01^ 

Since we are primarily interested in the kinetic behavior of f qq , 
we shall not specify the Q's and 7*s in this and higher orders 
in e (7^ is actually zero but higher-order 7‘s are not) but 
we will give the A n . We find 


dr P 00 

= - <£'(x,p) dss££'(x,p)>f 00 (2.61) 

The fact that this A^ commutes with V is' not made clear 
here, but this point will be covered later in this calculation. 

By using similar reasoning to higher orders in e , we 
have found 


3 f 


00 


37 


2n 


“ A n-00 


(2.62) 


where 


= <£'(x,p) J dss n £ n £'(x-ps,p)> (2.63) 


n 
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In Appendix B we show that 


n 


k-0 


k 


(2.64) 


where we have adapted the notation 

[£,<£ ] — [£,£ 3-^ j [£>[£,£ 31 = [<£,£■ 3g (2.65) 

so that [£,£'3^. stands for the k'th commutation of £ with 


In Appendix C we show that 


[£,£'] k = (-l) 1 / a ( k " 1 )r l (k odd) 

( 2 . 66 ) 

U,C'l k = (-ljVatk-Sjp^ {k 0ven) 

( 2 . 67 ) 

and, of course. 


[ £, <C ] q = <E 

(2.68) 

where and F 2 are defined by 


r, = [£,£'] - p-ln-n'-a'-o] —• 

1 dp 

(2.69) 

and 


r = [£,r 3 = p-L^'-p+p-fi'] 

dp 

(2.70) 


Equations (2.66) and (2.67) follow from the additional property 


U,r 2 ] = -r-L (2.71) 

In Appendices D and E we .prove that 
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(2.72) 


<.T 1 (x ) p)Z' (x'^y = -A(s)£ 


and 


<F p (x,p)£'(x,'p)> = A(s)(£ 2 -V 2 ) 


2 P,2 


(2.73) 


Thus,, we can write 

\ ■ (' 1 ) n£n ^I“ dsA (s) S n V 2 f 00 + 


l>- — n o U 


00 


k— 1^3 i • • • 
n 


- {i^7TW £n ‘ k I” aBA (3)3 n (£ 2 -V 2 )f 00 (2.74) 

k=2,4,... 0 


n 

Since £ and v commute, it is clear that each A commutes 


n 


AO 

with and therefore the compatibility conditions are satisfied. 

Prom Eqs. (2.25) and (2.29) we see that 


00 00 

3 £o v 3 £<x> n y s n A f 
3tT L 6 “ L 5 A n£( 


‘ p i — 1 ‘ pn 

2 n=0 2n 


(2.75) 


n=0 


We substitute Eq. (2.74) into (2.75) and in the terms containing 
the summations over k , we first interchange the order of 
summation, then substitute for the summation over n , the summa- 
tion over m = n-k to obtain 


bf 


^ = f°° <SsA (s)Y i (.e) n Z n ^ r s n V 2 £o + 

u n=0 

00 y / \ 00 


k=l , 3 , . . . m=0 

-j[ dsA(s)^ fes) k ^i^ £ m s m (£ 2 -.V 2 )f 

k=2,4, ... 


'0 


m=0 


0 

(2.76') 


9f> 



The summations can now be recognized as the expansions of the 
exponential and ' trigonometric functions. Thus 



f dsA(s)e~ e ^’ s [cos esV^-*-(l-cos es)£^- 
Jn 


sin 




(2.77) 


Along the restricted trajectory given by Eq. (2.38) we obtain to 
second order in e'n 


(p,t) + e£f 0 (p,T) = (eT dsA(s)e“ e< ^ s [cos esV^ + 

+ (l-cos es)£^ - sin es£]fQ(p,T) (2.76) 

Equation (2.78) gives a kinetic description of charged particles 3n 
a stationary random magnetic field . In the next section* the 
conditions under which this description is valid will be specified 
in more detail. 

3. THE RELAXATION TO ISOTROPY 


Equation (2.78) has the property of driving an arbitrary 
initial distribution function to isotropy. In fact, the final 
state is always the angular average in momentum of the initial 
state. To see this, we note that since £ and V commute, they 
must have common eigenfunctions. These eigenfunctions are the 
spherical harmonic functions with 0 , the polar angle relative 

to the mean magnetic field and 0 , the azimuthal angle about 

the mean field, as variables. We find with this choice 

$2y em - -MM-yijn (3.x) 


and 



-imY 


im 


We substitute an expansion of f^ given by 


(3.2) 
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(3.3) 


00 Z 


£q(Pj> t ) ^ ^ f 


Z-0 m=-i 

into Eq. (2.78)r and use the orthogonality of the spherical har- 
monics to obtain 


df» (t) 

157 “ 

=-(.sn)^J asA(s)e lm€s j^(i(^+l)-m 2 )cos ssm 2 -lm sin 6sJf^ m (T) 

(3.4) 


The effect of the integral coefficient on the right side can be 
made more transparent by introducing the spatial Fourier transform 
of A(s) . We imagine a sensor moving through the stationary 
magnetic field with constant velocity* U * measuring the 
magnitude of the components of the field orthogonal to U . The 
sensor sees the random part of the field varying with time. Thus* 
the Fourier transform can be written as a power spectrum which 
depends on the angular frequency* kU/X n * where k is the 
dimensionless wave number. Then 


AU) -i / 


dke 


iks 


.(f) 

v rr 

TJJUJ 


(3.5). 


and Eq. (3.4) can be written 


df 


- i(em-(eT,) 2 A' )f - 


(3.6) 


where 
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p a ( x 2 m ) *±(jr ( ra+1 )) 

“ “ 2 FjJOp + + 

P l(jr (m-l)) . 

+ |UU+l)-m(m-l)) (3.7) 


and 



+ 2 ( & ( 4+1 ) ■ -m (m+1 ) )(k +e ( m+1 )) + 2 ( ^ ( ' g+1 ) " m ( m “ 1 ) )(k+e(m-l ))] 

(3*8) 

3n which <P denotes the Cauchy principal value. Assuming a 
positive definite P ± * we see that is always positive and 

Aqo = 0 . Thus, all f^ m decay exponentially to zero with decay 

time (measured in units of the traversal time) 


'im 


(et}) A 


im 


(3.9) 


except for f qq which remains constant. We notice also that 
the f^ m oscillate with a frequency modified from the gyrofre- 
quency in the mean field by the factor (er}) 2 A' m . The differ- 

ence between the real gyrofrequency and the gyrofrequency in the 
mean field depends on the specific power spectrum which is in- 
serted into Eq. (3.8) . 

In order for the time scale extension technique which we 
have used to obtain this result to remain valid* the various time 
scales must remain distinct to preserve the ordering of terms in 
our power series expansions. Thus* the decay time must remain 
distinct from* and larger than* the traversal time. We must have 
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2 

a i5m ^ ( er l) « 1 j this condition is generally 

satisfied j ■ unless f^ is very anisotropic. In that case* we 
could have » 1 for large values of £ . Our kinetic 

theory is correct to second order in (erj) only if the more 
restrictive condition 

(sn)% m « i ' (3.10) 

is satisfied. 
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APPENDIX A 


We prove that 

<£'(x)£'(x')> = A(r)V 2 (Al) 

when 

x' - x - pr (A2) 

Jk Ay 

where r is the magnitude of the separation between x and x 
Prom the definition of we have 


<£'(x)£'(x')> 


- •ir-<&'(x).pp*ft'(x')>. 

dp ± 



CA3) 


but 




e lk^ mi n p k p m H in (f) 


where 


a 

r 


A y 

X - 


J 

X 


and where 


K M (r) = <P](x)^(x")> 
For. the isotropic random field 

R in (? > = A ^'> S &n- tB( - r)r tl r n 


(A4) 


(A5) 


(A6) 


where A and B are even functions of r . Notice that the 

second term., B(r)r*r , in Eq. (A6) does not contribute to 

. , . A y A /V 

Eq. (A4) when r - x - x is proportional to p . Therefore 


<n'(S)-pp.a '(5')> 1J - e i k ^ mJn p k p m A ( r ) 5 to 


= -(Aij-PiPj )A(r) . 


(A7) 
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Thus 


<£' (x)£' (x' )> = A(r )(Xc»n»^7r) = A (r )V 2 

^P 


where we have defined 


(A8) 


with 


V s = 

a P af 


n ij " 6 ij " 2 l*^j 


(A9) 


(A10) 
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APPENDIX B 


(a) We prove that 

*' £n -I 


where 


[£,£'l k = (-l) 1 / 2 ^ lc-1 V 1 (k odd) 
[£,£'] k = (_i)l/2(lc-2)r 2 ( k even) 

and 

U,Z'] Q = £' 


(Bl). 
(B2 ) 

(B3) 

(B4) 


The proof follows by noting that for n = 0 Eq. (Bl) is an 
obvious identity and for n = 1 and n = 2 Eq. (Bl) is also true 
if we remember the definitions of and T ^ 

r 1 = [£,£'] , r 2 = lz,r 1 ] (B5) 


We show that the (n-kL)'th term given by Eq. (Bl) is in agree- 
ment with the results of operating from the right on Eq. (Bl) 
with £ . Thus., any order can be generated from the n=0 and 

1 cases. 

Prom Eq. (Bl) 

£ ' £n+1 ■ I {nrfe- £tl ~ k [ £ . £ V < B6 > 

k=0 


Let n be an odd integer. Then 


3 ^ 



£'£ n+1 - £ n £'£ + 


+ 1 -fer^TEr £n ' k (-1 ) 1/s (k ' 1 )r i £ + 

k=l,3*** 

+f ^- k (-D 1/2(k - 2) r s £ 

k-2 , • * 


(B7) 


We define m = n + 1 and commute the and Tg which appear 

in Eq. (B7) to obtain 

£'£ m . A' - + 


>Z (-i)V2(k-i) £ m-i- k(£ri .r g ). + 

k=l,3, . . . 


i (-i) 1/2(lc ' 2) c m ' 1 ‘ k (£ r a + r i ) (B8) 


k=2j 4, . . . 

With k = k + 1 3 Eq. (B8) can be written 

m-1 


£'£ m = £ m £ 


- - £ m - i t£,£-] 1 +£ ■ f^4=j^re4 i £ m ' k u,£;] k + 


k — 1,3, • • • 


-I (-i) 1/2(K - 2) ^% + 

ic — 2 ^4^ . . • 
m-2 . 

. V (-1) (m-1) I P m-k rr . p'i 
+ Zj Tm-l-k)^T" * J k + 

k— 2 • • 


.(-!} 

i k_i i 

[m-1] 

>i 

(m~K 




k -~3j5j • * • 
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Now* letting fc = k and combining the summations over odd and 
even values of k * we find 


m-1 k 

£ ' £m - ♦ 1 Wrer £ "' k[c > £ \ + 

k-1 * 3» • • * 


m-2 


i n c - 

+ Z -fe^TWT £m ' k t £ ' £ ' ] 

1c=2*4* . . . 


k 


(BIO) 


or 


m , 

£ '£ m _ V (-1) iRl £ »-k [£)£ '] (Bll) 

k ^ 0 ( m -k)lk! 

which completes the proof for odd n . The proof for even n 
is similar. 
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(b) We prove that 


[£,£'] = Y ^ Lft , L- £ n ~ k £ ' £ k (B12 ) . 

k=o~( n " k ) Ik! 

for all n by noting .that for n = 0 and 1 , Eq. (B32) is 
obviously true,, and that given Eq. (BH2) for an arbitrary n , 
then the correct expression for the (n+l)’th term is generated 
by commuting Eq. (B12) with £ . Thus 


n , 

fr> (i'I . \ (“i) hi fr> p^"’^f' , p k l 

n+1 “ L (n-kjncr L<Lj,L Z z 1 

k-0 


(B13) 


Define m = n + 1 . Then Eq. (B13) is given by 


ir.~l 




k=0 

m 


*11 v 

Zj (rn-k } 1 (k-l ) l 


£ m ~ k £ £ k 


k=l 


(B14) 


'where the summation index in the second term on the right side 
'has been shifted by one. The summations on the right side of 
Eq. (B14) for k = 1 to m - 1 can be combined to give 


m-l 


[£,£'l m = £ m £' + (-l) m £'£ m +£ 




£ m - k £'£ k 


(B15) 


k=l 


which can be written 
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(bi6) 




m 


k=0 


(m-k)Ikl 


and which completes our proof. 
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APPENDIX C 


The T-operators are defined by the commutators ' 

[£,£.'] ■ ^ ' (Cl) 

t£ ' r i ] s T z (02) 

In addition, 

U, r 2 l = -r x (C3) 

To establish notation, we rewrite Eqs. (Cl), (C2) and (C3 ) in the 
following manner 

[£,£'].£ = [£,£'] = r 1 (c4) 

[£,£'] 2 b [£,[£,£']] = r £ ( C 5 ) 

U,z'] 3 = [£,[£,[£,£']]] = -1^ . (C 6 ) 

In this appendix, we prove that 

[£,£'l n = (-i) 1 / 2 ^ -1 )^ (n odd) (C7) 

and 

[£,£'] n = (-l) 1 // 2 ^ n- 2 ^r 2 (n even) (C 8 ) 

Notice that Eqs, (C? ) and (c8 ) reduce to Eqs. (c4 ) and ( 05 ) for 
n = 1 and 2 . To prove Eqs. (C7 ) and (c 8 ) , we assume them 
valid for the N*th term and prove that the (N+l)*th and 
(N+2)'th term are generated by the commutation rules of Eqs. (C2 ) 
and (C3). 
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Let n = N bs an odd integer. Assume 


U,£'] n = (-l) 1 / 2 ^ 1 V 1 • (C9) 

Then 

[£,£'] n+ 1 = (- l ) 1 / 2 ^ -1 ^,^] 

= (cio) 

Define M = N+L , an even integer. Prom Eq. (CIO) 

[£,£']„ - (-i) 1/2(M - 2) r 2 (on) 

in agreement with Eq. (C8 ) . 

Now 

U,£'J m+1 = (-i) 1 / 2 (M ' 2 ) [£,r 2 ] (012) 

Prom Eq. (C3 ) 

U,£') M+1 = (- i ) 1 / 2M r 1 ( 013 ) 

or 

t£,£'] N ^ - (.i)V2UN+2)-l) ri (C14) 

Thus, all odd n terms can be generated starting from Eq. (C4 ) . 
The proof of the even n terms is similar. 
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APPENDIX D 


We prove ' that 


<r 1 (x^)c'(*;p)> - -A(r)£ 


when 


x = x - pr 


The operator 37 can .be written in the form 


A -A 


' \ ^ 


r 1 .= p • 0'p-(3£') • 
Therefore, Eq. (Dl) is given by 


Sp 


-< p i( p i( x ) p r e i p j( x)) 

t 

which can be rewritten as 


, . r r x)) Pv e v OrfiL (* ' ) r*“> 
J J dp< 




~ a 

-Pi g~ p k e k£m| 




<L_ 


where r = x' - x . We assume the isotropic random field 


r i 3 ( ? ) - A ( i ’) 5 i J + B < r > r i r j 

Then Eq. (D5 ) reduces to 


dp. 


h 


-P- 4~ Pk s ktol (6 im A+ri:Bp iPm )P j “ -^i^jmA+^PjPm 5 


J a P.e 


'^PiPm^j = - p '' " ‘ “ ' * 3 ' 

Because of the skevj- symmetric € ^ m the terms containing 
do not contribute. Some further manipulation of Eq. (D7 ) 


(Dl) 

(D2) 

(IQ) 

(D4) 

(D5) 

✓ 

(D6) 

(D7) 

B(r) 
leads to 
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OPjg 


^3' 


+p i[( 6 jk“^j p k) + p k T^^ki^i f 


dp. 


*>i 


Again, because of the skew- symmetric e , only the first 
of Eq. (D7 ) contributes. Thus, 


< ri (i,p)£' (£;?» = Mr^c^gj 

, Xi 

= -A(r)£ 


(D8) 

term 


(D9) 
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APPENDIX E 


We prove that 

<r 2 (x,p)£'(x'p)> = A(r)(£ 2 -V 2 ) 

A ^ ^ A 

when x — x — pr by substituting in Eq. (El) the following 
expression for ? 2 

r 2 = (P-3')£ - £' 


Thus 


<r (x,p)£'(x,'p)> = £<(P*£"(x))£'(x'p)> - <£'(x,p)£'(x'p)> 


= £<(P’P (x))£ (x'p)> - A(r)V 2 


We must evaluate 


<(S.p'(j))c'(x;5» = -p m o m (S)Pie i1k p,;(x') V> 

a P;( 

€ i j k Pl mk ^ r ^ A 

a P . 

3 

± j , -i 

where r = x - x . With 


R mk (^) = A ( r ) 5 ™v + r 2 B(r)p m p. 


mk 


■nTk 


we find 

<(P*p'(x))£'(x'p)> ^ £A(r) 

Substituting Eq. (E6) into Eq. (E3), we obtain Eq. (El) 


(El) 


(E2) 


(E3) 


(E4) 


(E5) 


(E6) 
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CHAPTER III 


OMNIDIRECTIONAL INTENSITY AND FLUX EQUATIONS 
FOR ALL BUT "VERY LOW ENERGIES 



INTRODUCTION 


Most modulation phenomena occur in the interplanetary 
magnetic field in the vicinity € > 1 ; Thus, expansions in 

e , as presented in the previous chapter, are not well suited 
for studies of the modulation problem. Expansions in .1/e 
might be somewhat more useful; but, as we have already noted, 
such expansions have been carried out incorrectly until the 
present time. 

In this chapter, we present an alternate method for 
obtaining a useful modulation theory. We treat e ~ 1 right 
from the start carrying the full dependence on e in the trun- 
cated master equation. To reduce the complexity of the result- 
ing theory, we give up a full kinetic description of the cosmic 
ray distribution function by studying only the first two moments 
of a distribution function which is assumed to be very nearly 
isotropic. This approach is well suited to studies of the 
observed cosmic ray distribution function where anisotropies 
have been found to be very small. 
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1. THE OMNIDIRECTIONAL INTENSITY MB FLUX EQUATIONS 


We expand the distribution function about isotropy in terms 
of the spherical harmonic functions. The result is 

O 1 O /s. 4. t . 

p f = — I + jpjj: p * 0 + . . . (!•!) 

where I(t) is the omnidirectional intensity given by 

I(x) = P 2 Jdft p f(p,T) (1.2) 

j 

and 0(t) is the flux given by 

0(t) = p 2 J dfi p pf (p,t) (1.3) 

The integrations in Eqs. (1.2) and (1.3) are carried out over 

all directions on the unit momentum sphere. Notice that both I 
-1 

and 0 are differential functions of rigidity. The next term 
in the expansion in Eq. (l.l) measures the anisotropy in the 
cosmic ray pressure tensor, which we neglect. 

By substituting Eq. (l.l) into the truncated master equa- 
tion and using the orthogonality of the spherical harmonics on 
the unit momentum sphere, we find 

H = 0 (1.4) 

and 

|| - eO-0 = (er\) 2 J 6 % <3X<(n' •p)e"‘^ X (p.Q')>. 0 (T-X) (1.5) 

We must find the behavior of 0 (t) for large t and then, on 
adding weak spatial gradients in I , we can find the transport 
equation for I which is appropriate for the modulation problem 
The major simplification which is introduced by the study of 
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only the lowest moments of f is readily apparent from Eq. 

(1.5) • Equation (1.5) is in the form of a matrix exponential 
equation for 0 where the integral operator which operates on 

•» . /S 

' is no longer a differential operator in p . In addi- 
tion, because of the integration over SI which remains, only 
certain low-order moments of the operator \(U*p)e ^p*0 )/ 

contribute to the evolution of 0 ; the full complexity of this 

operator never enters. 

In the "interaction 11 representation of Eq. (1.5)* we have 

for 

S (?) = e~ ear 4(r) (1.6). 

the equation 

|| - a%r a^<(fi'-p)p-S'(x 4 )>-s(t-x) (1.7) 

where 

{^(x-r) = (x-r)^ - r 2 (x-r)cos eX + F^(x-r)sin e\ (1.8) 

and where 

? - [p + + < COS e x 6X - 1 ) ]-^ (1-9) 


To take advantage of the fact that Eq. (1.7) only depends on 
certain low-order moments of the operator K.{Sl'- p)p*fi (x-r))> , 
we introduce Legendre polynomial expansions of A(r) and B(r) 
from the correlation function. We let 


A(r) A i (X)P i (s) 


(1.10) 


4=0 


and 


B(r) 


( 1 . 11 ) 


4= 0 


where P.(z) is the 4‘th-order Legendre polynomial function 
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of z = (p*P) and where r is the magnitude of Eq. (1.9) given 
by - -s. , 


r = X 1 + (l-z 2 ) 


1 p 

■j|(eX) -cos ‘ eX 

i(eX) 2 


The coefficients^ A«(X) and B*(X) are given by 


\TL/2 


V X > - (HZ') 


^(X) = (^)£ B(r) V z > di 


After considerable -algebra, we find 


| . . (en) 2 ,3. J X dX^2J x oos cX-8X a K ( 1 )(2°£_J^) 2 ]p.S(T-X) 


(e-n) 


2 wJ 6X J x +J n cos 


- g ( s ° g .. 4 zl) . 1 k-s(t-x) + 

V (*x ) 2 ' 


(er ,) 2 T= fax J.. sin eX'-X^K^+K^ 2 )) • 


4 of cos eX-lVn 

2 v — — A 1 


lis^) n-s(t-x) 


where 
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and 


J ± - ^0 - I A a] (i- 16 ) 

Jjl = ^ ^ 4 2 + -o] (1*^7) 

K<11 Mst B 4 - f B 2 + B o] f 1 - 18 ) 

K (l1 + K (2) f ^ B 4 + | B 2 * B 0 J (1.19) 
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2. THE LIMIT, 6=0 


In the limit e = 0 (but (eu)^ finite) we find 

(!) 


K 


and 


15 "o * 


J 

°W 3 0 


rC 1 ) + r( 2 ) = b q 


(2a) 


j = i!£ a 

J -L 3 ^0 


( 2 . 2 ) 


In addition. 


Aq = A(X) and B n = B(X) 


0 


(2.3) 


Therefore 


j 

SS 

St 


= -2(er)) 2 f dX A(X) S(t-X) 
J 0 


(2.4) 


In this limit, 0(t) obeys the same equation as S{t) 

A straightforward linear time- scale extension of Eq. (2.4) 
yields 

J 


= _0 O 


(2.5) 


where 0 O is the leading term in a power series expansion of 
in powers of (erj) . We have, in addition, normalized all 
lengths to the parallel integral length 

X = 2 fd\ A(X) (2.6) 

P Jq 

Upon adding weak spatial gradients, we find 

|| + <*4 0 = o (2.7) 
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■©-U 



and 


3*0 

dr" 



1 

3 


a 

aVI 


( 2 . 8 ) 


In these equations, the gradients are measured in units of L , 
a macroscopic length which measures the size of the system con- 
taining the particles and 



( 2 - 9 ) 


For 

Eqs. 


J 

^0 

The result is 


l(*r) slowly varying in t , we can eliminate 
(2.6) and (2.7) when r > 1/a 

Sl/-V i i 1 

(vi n ) - i - 


0 + 
3t~ + 


a 


(en) 


— V^I 
2 v x 0 


from 


( 2 . 10 ) 


a familiar convection-diffusion transport equation for Iq , the 
leading term in a power series expansion of I in a . A new 
time scale has been introduced in Eq. (2.10) 


T = ctr (2.11) 

The velocity, V , is the average Alfven drift velocity due to 
curvature and gradients in the magnetic field averaged over an 
isotropic distribution of particles. This result is in agree- 
ment with the previous chapter if the expansion of Eq. (l.l) is 
substituted into the high-energy theorjr calculated there. 
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3. FINITE ENERGY (g ~ 1) 


Equation (l‘.15) can be written in the general form 

A 

^ dXK(ejX) *S (3*1) 

where K(e,X) is a tensor kernel which symbolizes the large 
collection of terms shown in Eq. (1.15). We have seen from 
Eq. (2.4) that when e 0 , K(e,X) takes on a very simple form 
with the attractive feature that it has a sharp cutoff in X 
due to the finite range of the correlation coefficient , A(X) . 

It is this property of the finite range of K(e,X) when e —> 0 
that allows us to apply the linear time scale analysis to obtain 
the high-energy kinetic and transport theories. When e ~ 1 , 
however, the kernel takes on a long tail which goes as l/X 
The mathematical techniques for handling Eq. (3*1) with such a 
kernel have not been developed previously. In the next two 
chapters, we develop the necessary techniques to determine the 
long-time behavior of Eq. (3«l) when the kernel is long-ranged. 

The long tail in K(e,X) is a general feature which comes 
from those particles which have pitch angles relative to the mean 
field close to 90° . The basic interaction between a charged 
particle and the random magnetic field lasts as long as the 
particle remains within a length X^ of its starting position 
at some arbitrary origin in time when the interaction is assumed 
to start. When e « 1 , the gyroradius of the particle is very 
large compared to X^ . Therefore, the particle’s trajectory is 
essentially straight through the correlated region; the distance 
between the particle’s position at time t and at time t=0 
grows linearly with time, and the length, X^ , is quickly 
achieved. On the other hand, when e ~ 1 , the gyroradius and 
Xp are equal in magnitude. Therefore, we must take into account 
the helical nature of the particle's trajectory in the mean field. 
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If the velocity of the particle parallel to the mean field is 
small, the particle may take a long time to get further than 
Xp from its starting point. The particle-field interaction 
can last very long and, as a result, can contribute an inordin- 
ately large amount to the interaction integral operator on the 
right side of Eq. (3.1). The result is K(X) « l/X for 
X » 1 .We present, here, a specific demonstration of how one 
of che coefficients involved develops the long tail for a Gaussian 
correlation function. By examining any one of the A^'s or 
B^'s one can convince oneself that the results shown here are 
a general feature of e ~ 1 even though A(r) and B(r)‘ may be 
short-ranged functions. 


- We ••••alculate Aq(X, e) when 


A(r) = e 


-r 


(3.2) 


We define the function 


C(eX) = 


1 - ^(eX) 2 - cos eX 


(3-3) 


Then 

r = X[(l+C(eX)) - SfeXJz 2 ] 1 / 2 (3.4), 

Therefore 


A 0 (^e) = e 


-X‘ 


: (i +IU\))J\ 


dz e 


-X 2 C(eX)z 2 




tt e ~^(l + C( € ^) ) — . , 

T e — Erf (VC ( eX J) 


x/TCixy 


(30) 


-X" 


Notice, for e = 0 , A Q (X, e) = e“ . However, for finibe e 
and for large eX , 
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(3.6) 


A Q (^e) 


— ^j(l - cos eX.) 



In the next two chapters, we study models of Eq. (3.1) to 
develop the mathematical techniques necessary to find the long- 
time behavior of the flux. In the first of these chapters, we 
consider a class of matrix exponential- equations in which Eq. 
(3.1) is contained. In the following chapter, we study a scalar 
model of Eq. (3*1) with a long-range kernel given by 

K(X) = T-irr (3-7) 

We will see that the behavior of functions governed by this 
kernel is rather surprising; in particular, the decay time is ' 
not analytic in the expansion parameter but behaves as 
(el Jin ej )“*** 
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CHAPTER IV 

APPROXIMATION METHODS 
FOR THE GOVERNING EQUATIONS 





INTRODUCTION 


We discuss here the uniform expansion of matrix exponen- 
tials in a rather general form. We succeed -in obtaining 
uniformizing formulae under the conditions that the generator 
of the time translations is a linear operator in Hilbert space. 
This means that we are not confining ourselves to a finite 
dimensional matrix system. The major restriction in the 
analysis is the assumption that the direct Taylor expansion 
of the unknown function has secular behavior which is of a 
polynomial nature. The polynomial has as its maximum degree, 
the order of the term in the direct perturbation expansion. 

This restriction is serious and effort is being made to remove 
it. It is, however, worth noting that with the definite poly- 
nomial assumption made very definite conclusions can be reached 
with practically no restrictions on the operators involved. 
Thus, a class of partial differential equations which is con- 
siderably wide is included. The type of equation included, 
in fact, -encompasses the class discussed by Akhiezer in his 
large monograph. 
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1. NONDEGENERATE PERTURBATION THEORY 


We write the equation to he solved in terms of a Schrodinger 
form as 



i || + AU = -SBU 

(1.1) 

We decompose the perturbation as 



B = OA + r 

(1.2) 

where we have chosen 

Tr (A • F) = 0 

(1.3) 

We readily obtain 

a = 

(1.4) 


Tr ( A • A ) 


Substitution of (1.2) into (l.l) shields 

i || + (1+05 )A • U = -sr . U (1.5) 


We redefine the parameters of relevance by introducing 

6 = (1+06 )T , (1 ’ 6) 

so that the original Schrodinger equation now reads 

i || + A • U = -er • U (1.7) 

We now proceed to perform an expansion of the following form 


b(t 0 ) = t 0 + ^ 0 (t 0 ) + ^S 0 (^ 0 ) + 

(1.8) 

Ti(t 0 ) = e^fTo) + e 2 ‘n 1 ('r 0 ) + ••• 

(1.9) 

T 2 ( T 0 ) = e 2 n 2 (t 0 ) 

(1.10) 
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We readily obtain for the derivative with respect to the para 
meter, s , 


d _ fds V 1 ^ , d , J2.f X d . « d V 

^ = Vd [3^0 + € ^1 ^ + € 1^1 ^ n 2 


(1-11) 


Similarly, we have 


Hr - 1 + s « 0 + 6 % + ••• 


Thus, the Schrodinger equation reads 


fd *d 2 .( '• d • 5 V 

+ 6 ^i + 6 V n i ^ + n 2 ^ ) + * * *_ 


. to = 


= “(l + € ^q + • • • ) (A + ef )U 


( 1 . 12 ) 


(1.13) 


We. notice that the expansion used (represented by Eqs. (1.8)- 
(1.10)) is sufficiently general to include both Lighthill 
stretching and a time scale. In particular, many time scales 
are included . 


In zeroth order, we find 

dU 


i ^ + .AU Q = o 


whose solution can be written as 

IAt, 


V*o> = e ’W 


(1.14) 


(1.15) 


A first-order equation can be written as 


dU, . dU n 

i + A .Di = -i? l7f 2 - i^-v 0 - r- u 0 


(l.l6) 


We can solve this equation as 



iat 0 a du n r T ° ~ 

- e O-TJ^) - ^ H- a 0 A • U Q + i J o dXr(X)*U 0 

(1.17) 
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where we have introduced the convenient abbreviation 


?(X) = e iAX -r • e~ iiU 


(l;l8) 


In order to study the secularity, It is useful to multiply the 
equation (1.17) by the Inverse of the leading-order operator as 
given by (1.15)* We then find 


iAT/ 


-IAt, 


dU, 


V^O^o) = e °*U 1 (- 1 )Uo 1 (^ 1 )e °- 3*2 u - 1 + i4 0 A + 

r T ° - 

+ i / dXr(X) ( 1 . 19 ) 


A simple manipulation yields 


-IAt 


IAt, 


T[ 0 = XI - iA-L + i a + 

p % 0 _ 

+ 1 / axr(-x) (1.20) 

n J o 


where we have introduced 


&» 0 ( T 1> -1,-* V 

A i - - 1 sir- 1 - W 


(1.21) 


¥e obtain 


-IAt 


lim e 
T 0~ >0 ° 


0 1 




iAr 


0 


= 0 


( 1 . 22 ) 


providing that the following condition is fulfilled 

A., = lim ■f'g— A + -i— f dXr(-X)'|' (1.23) 

1 t q »^1 *1 J 0 J 

We can make the calculation more explicit by using a complete 
set of eigenfunctions of the unperturbed operator. We let 

A|n> = a n |n> (l-24) 

We can then obtain matrix representations in terms of the basic 
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eigenvectors as 


A = Z a n^ n ^ 


n 


r =i <«i 


m,n 


We have then the following expression for F 

1AX 

1 ffl ,n |H!/ Nn ' e 

m 5 n 


r(-X) - e- lAX ^ r m)n |m> ^le 1 ' 


-i(a -a )X 

' m vn * \ 


-£ r m • ” n ' i*> <»i 

m,n 


-i(a -a )X 

> 111 n U> <nl 


We readily can split the sum as 

r( ~ X) < m l + E r m, n 

m m,n 

where the prime signifies m / n . We thus see that (l 
becomes 

Al % + irZ r mm^ m > < ra l} 

0 m 1 m 


If we now let 


and 


we can write 


*1 ~ T 0 


in = V* 


0 


A 1 =1^% + r mra )|m > <m l 


m 


Furthermore 3 

V A = I ^ a m + r mJ |ni> . <m l I a n l”> < n l 


n 


“ ■ y Z a m a tt>® > <ra| • f Z r .mn a m lm> <m| 
m . m 
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(1.25) 

(1.26) 

(1-27) 

(1.28) 

.23) 

(1.29) 

(1-30) 

(1-31) 

(1-32) 


(1-33) 



It is readily shown 


Tr(A 1 - A) = -y£ V“m - 0 
and we thus conclude m 

(1.34) 

O 

II 

O 

•> 

O 

II 

(1-35) 


The final expression for reduces to 

A 1 < ml C 1 ' 36 ) 

m 

This is one of the principal results of the paper. We notice 
that the first-order quantity is written as • 

a a r r r ° ~ . "1 .> 

Ul(T 0 ) = e u U.^) + i[ J o dXr(X) - t 0 A 1 Ju 0 (t 0 ) (1.37) 
Since we have 

[ a 1 , a 3 = 0 (1.38) 

we also conclude that 

. A x = A x (1.39) 

To obtain a final expression for the perturbation result, we 
notice some useful relations, in particular. 


'M 


r e 
mn 


1(a m~ a a )X 


|ra> <n| 


m^n 


r 0 _ v~"i / 

i aXT(X) = V 0 + X U 

u m,n \ 


/ 1( V a n> T 0 


- 1 


1 ( a m- a n> 


(1.40) 


|m> <n| (1.41) 


We can then write 


T 0 

J dXr(X) = +/(r 0 ) 


(1.42) 


where we have Introduced 


, i(a -a )X 

— V I'm n J * 1 

sS (\) =. l r m ( i(a -a 1 


|m> <n| 


m^n 


.m n J 


(1.43) 
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Thus, the first-order expression is given by 


W - 8 


U^) + il(t 0 ) • 0-«) 


Differentiation yields 


^ U 1^ T 0^ . iA,T o 


- 4(t 0 ) * A 1 - U 0 (t 0 ) (1.45) 


Multiplication by A^ yields 


j_iA X 

“ e ° ^ • V^) + IA-^Tq) * U 0 (t 0 ) (1.46) 


¥e thus obtain 


^l (T o) . A TT ,± , _ . iAT Or . A tt ri 


+ A 1 U 1^ T 0^ ~ 6 


A 1 U 1^ T 1^_ ( T 0^ A 1_ * u o^ t o^ 


The following calculation now indicates the behavior of the 
quantity ^ . We have, in fact, 

[J(t 0 ),aJ -X^.nC- T;"^) 0 ' 1 ) !, < n l’l“'> <ra 'l] 


(1*47) 


m,n m 


m 

/„ 1( V a n )T 0 

r r . — : 

, tnn m m V 1 ( a m - a n ) 


i)[|ni> <“>'IS nm -|m'> <n|6 m - m 


x „ x(a -a ) 

/ p «i n' 

r r i e 

mn nmV i(a -a 


T-^i)|itO<n|-^r mn r 


i(a -a )t a 
. x m n' 0 


a —a 
1 m n' 


^)jm><r 


= V r (r -r ) 

Z_, mn v nn mm' 


i ^ a m” a n^ T 0 


1 ' a m** a n' 


-)im> <nt 


(1.48) 


We then conclude 
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w5(T n ),A n ] =y r (r - r )(S— r 

v O'* lj i_j mn v nn mrrr\ x 


i(a -a )t 
v m n 7 


’(a -a T 
v m n 7 


m^n 


^■)l«v> <n| 


( 1 . 49 ) 
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2. SECOND-ORDER PERTURBATION THEORY 


The second-order equation is written 


i Sa 2 + A D 12 SU 1 SU 0 Su 0 

i ■*— + a • o 0 = -ie n 

OT 0 


J 2 - " l5 l ' li 'l 3 ^ " i '>2 3 ^ ' r ‘ U l ‘ V ' U 0 


SB 1 • • Su O 

- - 1 3?7 + Wo - ^2 XT - r - U l- V • CT 0 • 

1 2 ( 2 . 1 ) 


We have to use the following expression 

du. 


1 -&q * le 


iAt n 3 Ut (t,) r f'° ~ 1 

0 J o aXr ( X ) - T 0 A lJ i l U 0 < 2 - 2 > 


which was derived from first-order theory. Multiplication by r 
gives 


r-u^ = r-e 


IAt 


0 


i r x ° ~ 

1 (t 1 ) + ir dXF(X) •• U 0 - iT Q r* A x * U 0 (^ 0 ) 

(2.3) 


u 


'0 


The second-order equation can then be written as 


Su o iAx n dU,(x,) 

i ^ + A - U 2 = ~ ie 


dU,(x,) r T 0 - 

37 - —- + i J dxr(x)^u 0 - iT 

1 0 


0 A 1 U 0 + 


IAt, 


r z o ~ 

Ui(ti) - ir j d\r(x)u 0 + iT 0 r. a x • u 0 + 


- r- 


bu 


0 


Hr 


+ ^1 A 1 U 0 V ’ ~0 


(2.4) 


The solution of this equation can be written 
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IAtq dU-^T^) 


- T 0 p T 0 ~ 


" “ T O e ^ 


cju., i k 'r 1 ; n n w ~ 

+ J 0 dX 4 dX ' r W¥o( T o) * 


1 O Q P “0 ~~ lATp. j 

- 2 t 0 A 1 U 0 + 1 J dxr(x)e u U^) + 


tq 

- j r dxr(x) J^° dx'r(x').u 0 + J T T °cix(T 0 -x)r(x)A 1 u 0 (^ 0 ) 


“ iT) l A l U 0 cFrJ + lTl 0 A * U 0 
We now note some integral relationships, in particular. 


( 2 . 5 ) 


f °dx f T ° ax'r(x') = f °dx'r(x') = f T ° ax xr(x) (2.6) 

J 0 J o J o J o 


0 ~ 


'0 “X 

Tn - T, 


l Q t q I*q Iq 

- f ax f dx'r(x') - f dx(x n -x)r(x)= -x n / axr(x) ( 2 . 7 ) 
J 0 J\ Jq u uj 0 


We can therefore re-express the second-order quantity as 


U 2^ t 0^ T 0 e 


iAT Q dU n (r n ) 


°vv , r 

<5iT[ T 0 J q 


axrtxjA^tiQ) - | t|a®u 0 + 


+ i 


p 0 — iAT n a P 0 ~ P 0 , 

J dxr(x)e u^(t^) - j dxr(x) J ^ dx r(x ) • Uq + 


iT1 l A l U 0 * T '2^ + ir| 0 A ‘ U 0 


(2.8) 


We also note the following decomposition of the double integral 
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p U 'v p 0 — 1 p 'U ~ P'0 — 

/ dxr(x) dx'r(x') =i / axr(x) / ax'r(x') + 

Jn O* ^ dr> Jn 


0 ~ n'O ~ 
'0 J 0 


+ zf 0 °^j[V'[r(x'),r(x) 

= ^(a^q+sJKtq)) + 

+ I / °ax fV'[r(x'),r(x) 

^ J 0 J o L 

We can then write the second-order equation as 

au 0 


( 2 . 9 ) 


U 2 ( T 0 > = -^2 ^’"O ' iWo + - I T 0 A ! U 0 + 


- T 


O' 


A iAT 0 du-J^) 


iAr 


+ ■*-( A l T 0 +> ^ s ^ + 


37 : 


- | ^o+Vi^VV^q - I 1 °'i)pa>-'[r()‘').7(M 


(2.10) 


After some simplification, we obtain 

3 U 


W = “^2 377 


+ V ' U 0 ^ 1 A 1 U 0 " | T o[ A l' 5] U 0 


- T 0 e 


1 

2 


1AT 0r3u i (i. ) a 1 ~ iAx n 

37^ * iA l U l^l) J + 


p T o p\ 

*-v* «»— -» “ ■ 

/ dX/ dX 
>n J n 1 

r(\'),r(x)ju 
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As in first-order theory ^ in order to study the secular behavior 
we multiply by the inverse of the zero-order operator to obtain 


h V^o^o) = - 


&q v 0 + n 2 A 


i A - I _a 

1 n 2 i 2 n 2 l“i’"J 




T o 1A ' t or Su i ( V 


- iA i u i(-i). 


I i -iAT n 

Uq ( T l) e + 


• i-A T /-v v *i » “IATa -j -j 

t^ e 0 ^i> e -kh'*'* + 


r ° r x s 


/ dX/ dX 

J o J o 

r(x'),r(\) 


(2.12) 


We can express the manipulations needed in the following form 


" iA l T l _“ iAT 0 1 TT A \ TT - 1 /- n iAt 0 iA l T l 

e 6 t£ U 2 (t o) U 0 ( t 0> 8 8 


= -±a + — A - i — A. e 1111 1 

2 r\ 2 n 2 i n 2 

-iA. T- 


37: 


- i Vi("i)>‘o 1( "i )+ 


iA. t. 


+ 




i l ~ 1A i T i7 7 (a iA i T i l ~ iA i T i ~ iAx o 
2^ w * e - 2 n^ e 


J' °dx r (x'), r (*)]' 


IAtq iA-^r^ 
|e e 


( 2 . 13 ) 


The secular part is then obtained in the following fundamental 
expression 



r , in 

iA = lim \ -SJ A - i -i A 

Tq 00 L *2 


t o 

s 


n. 


-/C)U n (T, ) \ \ i \ t ~ 1 iA n T n 

Kafr- 1 - - <*1> + §M’> + 


*1 ~iA, T-, iAT A p T 0 p^- 

e 1 1 e °r ax/V 
^'*2 J 0 J 0 


r(x'),r(x) 


1At„ IA-iT, | 
e 0 e 1 1 (-(2.14) 


We will now perform several operations to simplify this result. 
We notice 


P T 0 pX 

-Vw 



T n 

r ® 

“ 

/ dX 6X 

b J o 

r(x'),r(\) 

L J 

= T o 


+ 2 ax 

J n 

^(X),A 
L -*• J 


+ 


+ / dxf^(X),3(X) 

Jn L J 


(2.15) 


so that we can rewrite the quantity Ag as 


f^Hn *l n 

iA c = lim A ~ i A- + 

2 T 0 -.J-i2 ’la 1 


T o "^Px 

-r~ e 
^2 


[( 5 ^ + ±t h ,J)} 1A ^ 


1 12 - -A e' iVl e'^ 0 

2 ”n 2 i n 2 


°cl^|5*(^),A 1 Je lATo e^ 1 ^- 2^f o °<3Xp(X),^(X) ,} (2.16) 


IATq iA n t 
le e 


0 r~ 


After some manipulations, we obtain 
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iA^ = lim 


- Ir 

— -> oo ^ *0 


, *1 


2 T 0 “» ^ ^ “ X 


A - i ~ A n + 
2 


T n -iA t, r/dU, (T n ) i \ t a , I iA n i 

“ t\^ e [G^ iA i u i( t i)) u o < t i> + ^ A i^l] e 


-iA, t.. -IAt^ p t 0 r~ i 

-% e s L ax r (x) ’ A i, 


1 IA-Tq iA.^-^ 
e e + 


n -iA t P T 0 r~ ~ 1 

5^ e 1 1 J o a^(x)J(x) 


iAi, iA.. t, 
e u e 1 1 


It is useful to recall the following relationship 


f~ I r->' / 1 ^%T%) X ,\ 

i(X),A n = ) r (r -r )(£-«-, <n| 

Lr 1J Z_j mn v mn mm'\ i(a -a n ) J x ' x 1 
m,n m n 


as well as the integral relation 


r X 0 r 7 1 V*' / o i(a m“ a n) T 0 , 

/ dXM(x),A.. = ) r (r -r )(- — - 1 

J 0 L i3 1. u m,n v nn mm'V ± 


m,n 


( a m“ a n ) 


T 0 \ 

„ -a 
m n 4 


TT- - IXa^O; 1 ^ 


The matrix representations for y$ and its derivative are 
useful 

J i x ) = ^ r m ,n( i(a -a ) ^)lm> <n| 

m^n 


-t- r — i 

- L 


i(a -a )X 

r ... e m n ' |m> <n| 


mn 


m,n 


(2.17) 


(2.18) 


<n| 

(2.19) 

very 


( 2 . 20 ) 


( 2 . 21 ) 
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The commutator of y$ with >/ can therefore be written 


W, ^ w s III li 


m,n m$n' 


X |jm> <n| , [m'> <n'| J 


= V" r r y 6 

Zl 1 mn m'n\ 
m,n 

m/n' 

- !■»'> <n|« n , mJ 


i(a -aJX 


m 


i(a -a ) 
v m rr 


n -iV 1 <V' a rf) x 


k ) 


[|m> <n'l« m , n 


|m> <Cn"| + 


y' r v r 

Zj mn Zj nn/ \ i (a-a j / 

nijti n' n^n' m n 

v - ^ V -1 / ^^ a m~ a n^ ,\ i(a ~a )X 

- I r m „I r m-m( S iT^- a „f -- ) e |m '>< n| 

m,n m' (2.22) 

mVm 


Some manipulations simplify this expression to 

~ I v~ 1 '’ r -1 / i(a -a )X 

- X I W- fifva - j ") 6 


|m> <m'[ + 


m^n m 

n^m' 


r-^'r-i / ^ a n** a m^ n \ i(a ,-a )X 

- ) > r T , (--■•■ . 7 ■ t~±)e m n |m'> <mi 

L L ™ m n\ i(a -a ) J \ 

m,n m' ' ( 2 . 23 ) 

m '^n 


We can combine the terms in the double sum in the following way 
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-^|m> <m'| + 


P«.®>] - I I {W.(^ 




m,n m 

n/m' 


>W 


- r r , , 

nm m n 


x(a ,-a )X x(a ,-a )X 
„ v m m' v m' n' ■ 
e -e 

i(a -a ) 
v n nr 


-)|m'> <mi j- 


(2.24) 


The secular part of the integral can now be extracted as 

... f «pM.5w] - £{r_r_ ifefl * ifeSj} 


m,n 


This expression can he rewritten as 


m n' 

(2.25) 


Sec f °dXrj(X),5(X)] = -2iT n y V T J m> ^1 (2.26) 

, L J 0/j mn nra (a -a ) 

m,n m n 

We are now ready to obtain a simple expression for the quantity 
Ag . Substituting (2.26), we obtain 


iA„ = lim A 

2 T 0 -> » W 2 >19 1 


+ 


Zo e - i Vir^ u i (:f i ) 
% 


. 6 ^ + [ v ,i } Iil ' 1 + 


x n -iA.T, rV / “^^ a m“ a n^ T 0\ iA -r 

+ 4 ° ZW r nn-WV^TI^-r-) m> < n|e + 

m,n 


m n j 


T 0 ” iA l T lV 

+ i n? e L r “ r 


jm) <mj 

mn~nm (a m -a n ) e 


m,n 


(2.27) 


We notice that this expression can be further manipulated as 



r % 

A = lim -{ -i ■— A - ~ A' 

2 l ^ r\o 1 


+ 


Tq “*> oo 


+ ^ + Vl (4 l) )r 0 i(i 1 ) + KVi]^ 1 ♦ 


T n ~iA n T 
- i e 


_ T-, r -1 ' / ^ a m a n^ T 0\ i 

1 l l r l m > <“l e 


1A 1 T 1 


n^n 


T q -iAj^T 


+ — ^ e 
n 2 


m 5 n 


!■»> < 


r r - pi , 

ran run (a -a„ ) 
m n 


iA n T 


l T lj. 


(2.28) 


It is important to express the commutation rules of A^ with >/ 
We obtain the following relations 


Ml - Z' I* « 


m,n 


(2.29) 


M = e ~ iAX [ A i^} 


iAX 


n' A _" 1 ^ a m“ a n^\ 

= -) r (r -r )(!=£, , — ) lm> <n| 

Li ran' nn mm'V i(a -a ; J 


m^n 


(2.30) 


We can thus write 


T *1 V nr • / -i ^ a m“ a n^ 1 \ 

K ,! f| - L r mn( r nn- r mm\— ~ - ) ^ < n 
m,n 


I (2.31) 


A simpler expression for A_ that results is given by 

d 
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lim 

T 0 ~* 





+ 


+ 


t -iA t r/ Su_ (t_ ) ' \ n j 

^ S [v ^ + + 


/T -F 


- i* < n ' 


m,n 


m iv 


iA l T l 
e + 


+ lH ~ iA l T l VV r < m 'i p iA l T l 

"Ho L r mn ran (a m >a n ) 


rn,n 


If we now consider the following two operators 


0 


1 


s 

V 

L 

m^n 


r (r -r ) 

mn v nn mm 7 

i(a -a ) 

K m n J 


|m> <nj 



m 3 n 


r r U> <mj 
mn ran (a -a ) 
' m n J 


we obtain the following relation 


[A r ° 2 l 


= 1 Vm'X TC^ ' O 1|m ' > <m ' L 


m 


nrij n 


m n 


-I Z W <ml6 ™-- U> <m ' l6 »'' 


m 3 n m , 


m n" L . 


which can be summarized as 

[A r 0 2 ] = 0 
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(2.32) 


(2.33) 


(2-34)- 


= 0 
(2.33) 


( 2 . 36 ) 



Similarly, we obtain 




' r (r -r ) 

mn ' nn mm' 


m 


m,n 


(aZ-a^) - 
v m n/ 


[|m'> <m'l , |m> <n| J 


■I Iw ( r "K , r ) ) [i”'> <»iw] 


m^n m 




' r (r -r ) 

m nn mm 7 


m 3 n 




[ r romi*^ > <n l ' l m > <n l r nn] (2 ' 37 > 


so that we have 


[■ 


4 i'°i 


/ 

-I 


m jH 


(r -r ) 2 

v nn mm 7 


|m> <n| 


(2.38) 


By writing 


* - A l u l^l> - | m > <n|U 0 (^ ) (2.39) 

1 v m n" 


m,n 


we obtain 


A p = lim |-i ^2 A ^ A, + r r . 

2 T -*«,! ^2 ^2 1 ^2^ mn M (V a -J 1 J 


m n' 


We are now prepared to determine the scales and we introduce them 
through 
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( 2 . 41 ) 


n 2 - t 0 

n 0 = ^o x 0 

^1 = 7 i T o 

We can thus write 


-i7 0 A - 


¥1 


r. |m)> <(ml 
+ Zj r mn r nm (s^-a j 
m,n 


This expression 'is more compactly written as 

^2 = ” i7 0 A “ 7 1 A 1 °2 

Multiplication and the tracing operation yields 
Ag ♦ A = -i7 Q A « A - 7-j^ * A + Og • A 
Tr(Ag*A) = -i7Q Tr A • A + Tr(Og’A) = 0 
We therefore obtain 


(2.42) 

(2.43) 


(2.44) 


(2.45) 


(2.46) 

(2.47) 


Tr(Og.A) 
7 0 = " i Tr(A-A) 


(2.48) 


Multiplication and tracing with Og yields 


0 


2 


A = 



r 

mn nm 


fii,n 



ln'> <n'| 



m,n 


r r 

mn nm 


a m |m> <inl 


(2.49) 
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(2.50) 


Tr 


(o. -a ) = V r r (- ^ 

v 2 1 Li ran nm\a m -a n / 


nun 


Thus, we can write 


7 0 = -i 


r r V 

mn nmV- w 
t, i. v m n' 


(— J sJ) 

rt V a m _a «/ 


I 


m 


a 


m 


Multiplication and tracing with A 2 yields 


A 2 * A 1 i7 0 A * A 1 " 7 1 A 1 ' A 1 + °2 * A 1 


Tr(A 2 *A 1 ) = 0 = -7 1 TrfA^.A^ -i- Tr(0 2 ‘A 1 ) 


and thus we obtain 


Tr(0 2 .A 1 ) 

7 1 “ Tr^-Ag) 

Multiplication and tracing with A^ yields 

A-, • A- = ) T ) T , ,|m> <m||m'> <m1 

1 1 L, m/j m m 1 x N 

m m' 

= 7 T r |m><m| 

/ . mm mm ' N 1 


m 


Tr <W=I r , 


,2 

mm 


m 


Since we also have 


(2.51) 

( 2 . 52 ) 

(2.53) 

(2.54) 

(2.55) 

( 2 . 56 ) 
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o * a = V r r .‘ >-<4 \ r , m s /m ' 

U 2 ^1 Zj run nm (a-a) m m ' ' 
m.n m n m' 


m,n 

-i' 

m,n 


r r p ■ l m > . < m l 

ran nm in (a m ~a n ) 


(2.57) 


^(cyV » Y. 


lUjll 


r r r 

mn nm mm 

(a -a ) 
v m n y 


(2.58) 


we obtain 


I' 


_ ro,h 


r r r 

ran nm ram 

(a -a 5 
v m n J 


1 


m 


mm 


We can then write 


(2.59) 


Tr(Og*A) • Tr(0 2 ‘A 1 ) 

2 = ° 2 “ fr(A^A]r A “ Tr ( •0 1 ) A 1 


(2.60) 


0 


2 


V' T r < m [ 

■ L mn nm (a ra -a > 
m,n 


( 2 . 61 ) 


This is a rather compact equation for the second-order unifonn- 
ization. We can now verify a number of the properties of the 
second-order uniformization. In fact, we have 
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Yr ^ n l 

„ A mn i(a -a 


m n 


r ,A i_ 


U UV a rJ ^ 

111,11 m 


]T r m 'm' ImXnl „ lm'Xm'1 


m,n m 


I Vm{ |m><ra '' Vn-I ra Xn l 6 m'm] 


i (a“-a n ) [ r nn |m > <nl-l">>ln> r J 
' ni n' L - 1 


TT^O ( r m- <“l 

m n' 


(2.62) 


Since the first-order equation can be written as 




7 p 

+ A-JJ^t^) = +i_Q] Ijs^rj < n l^i_ u o( T l^ 


ra n j 


( 2 . 63 ) 


we have 


; iAlTl V^) +/ Tl dT'e iAl(Tl ’ Tl) 




n> <nQ,A x U 0 (t') 


iA n r n . pT- iA,Xr/v— 1' F 

■■ '■£ “* *[(1tv€ 

u m,n 


\ ~ -iA n X 

•|m> <nj J, iA^Je “ U 0 ( ^ 


( 2 . 64 ) 


After some manipulation, we can write the result as 
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* 1A-, ir n . iA. t. \r- >' F — iA. ' 

V T l> -e 11 0 1 (T 2 ) + e X1 X (a-4) l m >< n l e 

m,n m n 


-Z'iw^< nl ^ 


m,n 


If we now Introduce the useful notations 


A l |m> = I <“» “ r ,J m> 


m J 


A ll ra> " r nJ™> 


we can write 


* iA-r,. iA.T. . ri'/ r \ 

U 1 (t 0 )= 8 e U 1 ^ t 2^ 


m,n 


x(r -r )t, x(a -a 
„ v ram nn ; 1 „ ' m n 
e e 


mn 


m^n 


a -a_. 
m n 


|m> <nju 0 (T: 0 ) 


We can then give compact results which are 


iat^ iAt-, .a vy r„„ Y e 1 (a m ~a n )T 0 ^i(r mm -r n )T 1 ^^ 




m^n 


m n' 


v iAT n iA n T 7 iA t , 

D q (t 0 ) - e 0 e 1 1 e 2 2 " 


VV 


A slight rearrangement also yields 


1 "(A* + 

(2.65) 

( 2 . 66 ) 

( 2 . 67 ) 

°im><n|U 0 (r 0 ) 

(2.68) 

|m> <n|u 0 (T 0 ) 

(2.69) 

( 2 . 70 ) 
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V* 1> - el4lTl u l^ a ) + XG%-)[e l(rmm " rnn)Tl - l]|nv> < n |a 0 (x 0 ) 

^ x m n /L J 

(2.71) 


m^n 


We note that since we have 


A 2 * °2 * - Vl 


with 


■r 


r r 

mn nm (a -a ) 
v m n' 


( 2 . 72 ) 


(2.73) 


we can also obtain a matrix representation for 0^ 
given by 


m/n' 


i \ V r- \ m 'y <( m i i \ 

- 2j r m'n' r n'm' (a m ,-a n ,j ^ 


This is 


If we then set 


we find 


■I 

n 

n^m 


r r |nl> 

mn nm \ 

n 1 m V 


=X?w lm> 


n^m 


°2 lm> = °iJ m > 


0, 


m 


-l 


r r 

m n nm 

(a -a ) 
n m xi' 


n/m 

We can express then as 

Agt 111 ^ - (0 m - ^o a m ~ ^l^rnm^ 1 ^ 


(2.74) 


(2-75) 


( 2 . 76 ) 


(2.77) 
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which can be summarized a3 


4>> 


= 6 |m)> 
m' 


( 2 . 78 ) 


°m x ^O a m ~ ^l^mm 


(2.79) 


We can finally rearrange the first-order result which can be 
given as 


>? 0 ) iATQp dU^(T^) 

— + A l U l^o) = 8 1 


_ i OT^ ~ + A l U l( T l) ]“ i [^(' r o^ A l_ * U 0^V 


= is 


x ^ a m -a n' 


|m> <n|). 


A lj e U 0 U 0 ) + 


pi' / i ^ a m“ a n^ T 0 ,\ A 

- i ) r (r -r m> <n|u n (x n ) 

mn v nn mnr\ 1 ( a m “ a n ) / 0 V 0 J 


= ±y .,V- , (r -r )e i(am " an)T °|Bi>^|o,,(^ + 

/, i(a -a ) x nn mm/ . 1 ' v 1 O' 0' 

J~L ' m n J • 

m*n 


, , i(a -a^)^ n 

Z A, m n ' O -, \ ■ » 

r (r -r )(£— — — — — )|m> <n|u n (T n ) 
mn v nn mm^ i CaT-a^ ) ' x 1 0 V 0 y 

x v m n' 


- 1 1 r m n( r „n- r mm ) r|#7 U 0<*0> 


r/<-*' r 


- HvA nw |ra 

m,n 


> <ni^A 1 J.U 0 (T 0 ) (2.80) 


which can be written compactly as 




+ A 1 U 1^ T 0^ 


x ' a mT a n J 


n> <n|^A 1 U q (t 0 ) 


(2.81) 
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Utilizing Eq. (2.32)* we now notice that alternative expressions 
for Ag can be given 


A, 


f n 0 

= lim •{ -i — A - — =■ A n -l- 
r Q -«* l ^ ^P 1 


2 o e 1 Vi 
^2 


/ Su-.CtO * Y , ^ 

X 1 ^ * A l U l^l)>0 ^1^ + 


m 3 n 


-Z'^nT^ !■» < n lK Vl +%T V- 1 ^ 8 - 


iA-, t- 


lim j-i 

oo 


Tq^CO 


- i 


* % 
^2 


m,n 
1A-, T n 


m rr 


* - "i \ ^ lTl * 

r«' r r* -j -i(r -r ) T , n 

1 iT 5 ^tI^<"Ia} ” nn l J- 

+ ^ , r_r_ _ ( 2 . 82 ) 


r 


n L mn ran (a -a n )J 
m*n 


We see thus that the second-border theory can be made fully 
explicit . Operators that uniformize the expansion are computed 
explicitly and the major assumption* we repeat* has been the 
linearity of the perturbation expansion in first order and of 
the quadratic behavior in the second order. 
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3. DEGENERATE PERTURBATION THEORY 


We now consider degenerate perturbation theory by writing 
the Schrodinger equation as 

i H + (A+eB)V = 0 (3.1) 

If we consider the nevj operator 

V = e iAT U (3-2) 

we thus obtain an interaction representation 



= -eBU 


(3.3) 


with the operator 

B = e'" lAT .B.e :LA ' r (3*4) 

We now introduce matrix representations as follows 

Alm,a> = a m |m,a> (3-5) 

The symbol a denotes the degenerate eigenstates associated 
with each eigenvalue a n 


B 


B 


I 

mn 

af3 

V 

L 


B ma,npl ma > < n ' P l 


B ma,nP e 


-I (a -a 
' m n 


)t 

|m«> <nPl 


m } n 

a,p 


(3.6) 

(3-7) 


The relevant time-dependent operator given by Eq. (3-7) can be 
rewritten as 


-X |ma> <mP1 

m 

00 


K — i _ £ *3. ) TT 

+ e m ” l ma > < nP l (3.8) 

00 
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Itfe now consider a linear time scale extension which is given by 


cHJ 0 

1 = 0 » 

. °J.O 




The zero in first-order equations yields 


dU n dU, 

1 + 1 ^ = " BU 


0 


or 




du. 


'1 v 0 
1 ~ - 1 - Ba o 


We can integrate to yield 


v A 8U n n ’U ~ v 

iU i( T o) “ iU l( T l) “ iT 0 "StT “ ^(^)U 0 (ti) (3-12) 


0 — 


(3.9) 


(3-10) 


(3.H) 


To investigate the secularity^ we multiply by the inverse of the 
zero-order operator obtaining 

M T 0> U 0 ( t 0> = U l^l) U 0. ^ T l^ " T 0 <5r^ U 0 ^ + 1 d ^ B ( X ) 

(3-13) 


We now define 


A- = -i ttt 

1 OT 






(3- 14) 


We thus obtain 
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_ U 1^ T 1^ U 0 " iT 0 A l + + 


m 

aB 


+ i ) B. 

m^n 

oB 


; (*— 

ma,nBV -i 


- i ^®ta- a n)‘ r O 


^ a m“ a 7> 


—* )|ma> <nBl (3 .15) 


¥e let 


A 1 “ I B m a,mP |ma > <mPl 


( 3 . 16 ) 


m 

a, 3 

so that the first-order operator is given by 


D l(*o) ■ W - Z B ma,nP 

m/n 
aB 

We note the commutativity 

[A 1 ,A] = Y,\a,mP [|ma> <mB| ,A] = 0 (3-l8) 

m 

aB 

The following matrix representations are used extensive ly 


( a m- a n> T < 

a m“ a n^ 


- i mct> ^BIUqC^) (3.17) 


A 1 

-I 4 ? * 


i 

B ma,rae |ma> <raPl 

(3-19) 


m 


aB 



A 


> 

S 

II 

a 

m 

\| m 0 /* \m^l 

(3-20) 

- 

m 



a 


We notice 








[A^A m ] 


0 

(3-21) 

We thus have 

a compact 

version 

Of 

(3.. 17) as 



83 



VV = + 1 / ^ < 3 AB(X) - T^jupC^) (3.22) 


We now define 


>/('r°) = dXB(X) - t q A ;l 


(3-23) 


The matrix representation is given bjr 


-i(a m ~ a n ) ,T o 

^( T 0> = *1 B m a, ng ( e a -a <nPl (3 ' 24) 

JrL m n 

m^n 

af 3 


The first-order operator is then given by 

W - 

The time derivative of Eq. (3«23) is 

= B(t q ) - A 1 

We can therefore write 

+ ±P\T; n ) 

'1 UL 1 


du x du^) g 
T€T = ^ + W T o> ff: 


or 


5u i _ StJ A) /r u TT , 

FrT " FF; ^( t 0 ) A 1 U 0^ t 1^ 


1 w ‘l 

resorting to the matrix representation 

-i(a m “ a n ) T o 


- X 1 v, 


m^n 

ap 


n{3 


|mq> <n£| 


(3.25) 


(3.26) 


(3.27) 


(3.28) 


(3 <-29; 


whenever necessary. 

We now consider second-order theory. The basic equation is 
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dU n dlL bv 0 ~ 

1 <zq + 1 ^7 + 1 ^ - - BU i 


(3-30) 


Using the lower-order results, we obtain 

Su 0 au - 
" - 1 3^ - 1 - BU 1 

3u n 3 d (t, ) „ , r \/ v <j \ 

= -1^-1 +i/(T 0 )A 1 u 0 (T 1 )4^(T 0 )+A 1 ^u 1 (? 1 )+iy)(T 0 )u 0 (-t 1 )) 


- -1 


. SD 0 s 3U l^l) 


,Jk 


+ iyf('r 0 MiU 0 ("i)-^^o) IJ 0^l)“ i ^ (^K) )U 0 ( ^ + 


-A l U 1 (' r l) ” iA iVJ(T 0 )U 0 ('r 1 ) 

• +A 1 U 1 (T 1 ))+ik/(T 0 ) ; A 1 Ju 0 (T 1 )-/(T 0 )U 0 (^ 1 ) + 
-^(Tq^Tq^qC^) ( 3-31) 


du Q / * 0 ^) 

1 77 = \ 1 


m 2 V o Tl 


whose integral can be expressed .as 

dU, 


* dU , dU, (t, ) j \ 

1U 2 ^ t 0^ = 1U 2^ T 1^ “ T 0 IL ck“ “ T 0V 1 + A l U lKv 

T 0 

+ i J dX[>/(X),A 1 jU 0 ( : r 1 ) (tqJUqC^) + 


- i J o °dX>}(X)vf(X)U 0 (T 1 ) 


(3-32) 


We shall need some integral representations. In particular, 

Tq » Xq 

f ax^{x)*f(x) -if dX[>/(X)/(X)+Kf(x)^(x)] + 

Jq d Jq 

■$r° 6 XV(X)J(X)) 

r • 

•|/(t 0 ) +|J ° ax[j/(x),/(x)] (3.33) 



i , dU n / Su i ( t, ) J \ 

iu a (T 0 ) = iu 2 (t 1 ) - V 3^2 . i. 1 - + a 1 u 1 (t :l )J+ 


+ i J °dX[^(X.),A 1 ]U 0 (^ 1 ) -J(t 0 )V 0 (1 1 ) + 


- 1 |ii 2 (T 0 )u 0 (^) - 1 1 r°axU(x) ;l j(x)]u 0 (^i) (3.34) 

To study the secularity, we use the trick of multiplication by 
the zero-order operator. We obtain 


f ^ \ tt r ± \tt“1 f - \ _ U 0 (T 1 0*1 ^ + 


. -1,-i . , -J . -1/-1 V ( 

W 2 (t 0 )V< t 0> - V^C K> - T 0 cir) 


. (. 3d i(" t P 


+ I'CA i 


0\ OT- 




+ f d\>^(X),A 1 J + ±J(t 0 ) 


- ^ (*r 0 ) - | f °(3x]j(X),J(X)J 


(3-35) 


The commutator of j/ with A 1 can be computed explicitly. The 
calculation yields 




B ma,rip(' 


-i(a -a )X 
e m n -1 


m^n 
a£ 


ma s nPv a m “ a n 


B_ 


m n 


Z Z B m^ nP B m 'p, m - vV a m -a 

m^n m' 
a£ p,v 

- |m'|i><nP|a m , m 6 J = 


m'P.,m'VL 
m' 

M-V 

“*-( a nT a J X 


m.c£KnP | ,|m' [jCKm ' v 


m n- _ x 


|ma><m'v|5 I)1 , n 6 (1 p + 
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1 X B ma,np(' 


-i(a -a )X 
e m n -1 


m/n 

ap 


a -a 
m n 



B nP,nv ima> <nv| + 


- B I m 

m|jL,ma 


v> <nP|} 


1 Yj B ra“, nP B nP, nv(' a m - a n 


- 1 <V a n) X 


m/n 

ap 

v 


— ^ |ma)> <nv| + 


1 Ij B na,mpC' 


-i(V a m) X 


-1 


a -a j B nv,na l nV > < mg l 

n m / v> 


m^n 
ap 
v 

-i(a -a )X 

A m n -1 

1 Zj B ma,np B nP,nvV a_-a” 
m^n 
aP 
v 


lma> <nv| 


+ 1 X B nP,inoPnv, 


m/n 

ap 

v 


c , 1 ( a m' a n) X 

nPV 1 ) I nv > <ma| 

KV m n 


The integral of the commutator is also of important 




Z. B ma,nP B nP,nvV 77 _ a \ 2 

m^n v m n ; 

aP 
v 


-i (a m -a„ )r n 
_ ' m n' o - 

— |mct)> <(nv| 


X B nv , nP B nP , ma(j 


m^n ^ a m _a n ^ 


i(a -a )t 0 
* m n' u i 

— j nv> <ma| 


ap 

v 


(3-36) 


+ 


+ 
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, T V B B l ma > < n ^l . 

0 ^ inc^nP^n^nv (a -a ) + 

m/n m n 

op 

V 

_ , T Vr, R l n V> < ma l 

0 4 ^nv^n^np^ma “"(a^-a n ) 


m/n 

ap 

v 


(3-37) 


The commutator of vO with the time derivative is also computed 
explicitly. We first obtain 


[i(x),j(x)] - -y y 

L 1 L 1 


-i(a -a U 


*' i ^ a m‘“ a n )X / e 

*0 Oi T5 1 ° 

mo^nP pp.,qv\ a„-a 


P q- 


m/n p/q 
op p.v 

*[|ma> <nP| „ | PM> <qvj] 


P q 


=-n 

m/n p^q 
a£ M-v 


B ma,nP B pp.,qv e 


-CV-We' 1 ^' 3 ^.. 


q p _a q 


[ |nia> <qv| 


Wfti - \W> <^i 5 qm 5 va 


m^n q/n 
ap v 


V „ aq)X -l\„A/„..i . 

Z_. B ma^nP B nP,qv e V a -a *yl ma Xqv|-f 

n- ^ 


’T’ r 1 -i(a -a )Xx ” i ^ a p -a m^ 1 \ 

+ X 1 B m a,npV,ma s (~ "a'-a, ~X XnPl 


m^n p^m 

aP p. 


V T ' 1 ( a rtr a n ) Ve" 1 ^ n ' ap 1\, N , 

‘ 1 l B ma,nP B nP,m e 


m/n p/n 
ap |i 
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+ X X ^nPjmoPpix^nP 5 


-i(a -a )\/ ( a p” a n^ X 
v n m' / e * 


m/|X p/n 

c£ m. 


V a n 


|pH> <ma| 


(3.38) 


Some manipulations are needed in order to write this expression as 


X X B ma,np B n0, 


m/n p/n 
a£ |x 


z -Kv-V*- -KvSiK 
«,(- ^ >*> <H*I + 


m^n p^n 
ap • |i 


/ -^VV 51 '^VVN, 

B W,,ntPnfi,ma(- 5^ ■)»> <ma| 


m/n p^n 
ap jx 


B B ( 

ma,np nP^piJ-V a p" a n 


-i( a m“ a n) X - 


■) |ma> <p[JL | + 


m/n p/n 

a£ p. 


-i(a -a )X 
v n m' 


-I I V,neV,mc(--- a -=F ; - ) l»> <““l + 


tti^h p^h 
op p/m 
M- 


p n 


•‘(vV^ 


+ 1'1 \a,n^,Jr^FT—) < p ^ + 


P n 


m^n p/n 
ap p/jn 


B piA,n£ B nP,nic\ a p -a n J ^ + 


lma> <m|xj 

ina^nB^nPjmix a m *' a n 


+ > B, 


|mix> <mal 

mix , nP n0 , ma a -a 


(3-3.9) 
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We also will need 


<3Xy(X)J(X)] » 





in n / 0 ,\ 


1 ( a m- a n) T 0 


m n y N p rv 


'T 1 'T 1 /p^^W^P^O ,\ 

+ X L L B ma,nP B nP , PH\( a ~a ) (a -a j) 
m.?n pfn F y 

ap p^m 

M- 

i(a -a ) 

V* V /<= v m p' 0 ,\ 

Lj Zj B pp.,nP^nP,mcAja“-a p ) (a -a // 
m/n pfn v v 

ap p^m 

M' 


+ T 0 B ma, nf^nP^mM- 
m^n 
apix 


|ma)> ^mM-l 
a m -a n 


+ 


V t-. -r, l mpC> <mct| 

+ iqX B m[Ji.,nP B nP*ma a -a 

m/n m n 

Qpp 


(3-40) 


We note that we can write 


SU 0^l) ~l,-> s iA l T l ^0^ n -l/A v ~ iAT l 
U 0 ( T l) “ e cRI U 0 


(3.41) 


We now define through the equation 
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i 


(3.-42) 


a °o ( A D -i,i j _ 
U 0 ' T 2' 


-A, 


This can also be written as 


s «cA> „-i,- , 

3^ U 0 (T 1 ) - e V^ 2 


1 l i “iA-i "^T 

(iA 0 )e 


(3-43) 


We thus obtain 


i T i. A Su i^i) 


* *V * * +lT ' 


T 2 

•u’ 1 ^) + i^(T 0 ) - (t 0 ) + 

T • 

“ 1/ °dM^(>.)^(x)] (3.44) 


The secular part of this expression can be extracted explicitly. 
We have 


a _i a 1A-.T, ~iA n T 

Sec ( t q)^o ^0 ) = ~^ T 0 e Ag® 


. 3u iW> 


l T 'l f 

W* 3^ 


+ 


A 1 U iA>X 1( A • lT ol B ma,nP B n(3, 


|ma> <nvl 


rn/n 
exp 
v 

- iT A ) By,., „ftB. 


nv 


lny> <jna| 


(a -a ) ~ 1T 0 A nv,nP nP,ma ( a _ a ) 

v fll XX * — t K m -n J 


m fn 
op 
V 


m n- 


B _ „ qB.. 


jmcx> <iW|x | 


2 A ma,nP nP,mn a m ~ a n 
m/n 
aPu 

T 0 V P P lm(i> < ma i 

2 A m{i,nP nP,ma a _ a 
m7n m n 


(3.45) 


opp. 
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The time dependence of A^ that is obtained in (3*43) is then 
obtained as 


iA,T 1 ' -1A, t-, 

e’ 1 x A 2 e x 1 


f aojdj) 

V 1 vq — 




+ 


^ ^ im«> <nv| 

* Li B mo^nP B n{Xnv (a -a ) 
m^h m n 

o£ 

v 


' X B nv,n£ B nf3,ma 
m^h 
ap 
v 


1 nv> <maj ± V' |ma> <m\ 

(a -a ) + 2 A ma 5 nW } ran ( a nf a n ' 

v m n' m/n 
ctPii 


i V. -r, |mp> <mal 

+ 2 Zj mp-^nP nP,ma ( a -a ) 
rn^n m n ' 

o0|x 


(3-46) 


We now introduce a new basis as 

!ma>* =^C^ y |m7> (3-47) 

7 

where 

A 1 |ma>* = 6^1 ma>* (3*48) 


The degeneracy is removed in the new basis by the fir.st-order 
perturbation* We note the following relation 

A}ma>* = a m ^C^ 7 |nry> = a m |ma>* (3-49) 

7 

The orthogonality of the new basis is expressed by 
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(3.50) 


*<nP|ma*)> = ^ Cp C^<mixSm7> 


7V 

\ ^ p^fl p^ 

- L Way 


In view of completeness, we can write 


2 _ ptn pin 

6 op "Z^cry S-P 

7 


(3-51) 


whence 


I 

cry 


+ + + 
c m c m c m _ c m 


'\ia ''ccy ^P “ v |xP 


Therefore, we also have 


n r> I 

5 = > C m c m 

\iy ix a ary 


a 


We can then write 


(3.52) 


(3.53) 




a 


a y 


“ZV-y |m7> = |mli> 


(3.54) 


The new basis can be written as 


;m|x> ! ma > 


a 


The adjunct vector is then computed by 


(3-55) 


<m\ * ^ 


tfL < ma l 


* 


a 


=Xc* <mai 


a 
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(3.56) 



so that we have 


<m|x| 


( 3 - 57 ) 


a 


The time-dependent A^ can then be written as 


iA 1 T, -iA,T-. 
e 1 e 1 1 


■ g% ( ti) + vi^i)X i( 4 i ) + 


S' B B i ma >< nv l _ V B B l mV X ma l . 

L ma,nP nP,nv (a -a ) A nv,nP nP,ma t _ a \ 
mfn m a m^n v m n' 

a|3 a ^ 

v 


V 

i 


X B ma,i 


C ra V 


TO i mv > # *<^ + 


2 A Z_j ma,np np 
m^n vy 
apjx 


m n' 


CmIc ? 1 


+ t L L B mM- J nP B nP,ma l mV > <*?l (3.58) 

m^h vy m n 

aPli 


We now assume that the starred basis is made by nondegenerate 
eigenfunctions of A^ . We can then let 


, ad =Yn n i ma > < nv l TT (i N 

dT 1 1 U 1 A B ma,n3 B n£3,nv (a -a ) U 0' T l' 
mfn m n 

ap 

v 


Y n 13 I nv > -* 

U B nv,nP B nP,ma ( a -a ) U o( T l^ + 
' ' m n' 


+ „ a B. 

m^n 
«P 
v 
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^ ^ c m c m 

- II ZA^nfsBnP,^ |mv>* *<xay\v 0 (\) 


m^h v^/7 

aPp. 


m n J 


o m+ c“ 


^ 1 r 1 ^ (j 

2 2j Zj B m(j. , nP B nP ,ma ( a - a ) ^ m ^l u o^ T l^ C3*59) 


m/n v/7 
apy. 


ra rx 


where the quantity A g is given by 


p m + pm 




-f- 


m^n v 
apia 


o^c 11 


+ 1 1 E B mlSne E nP,ma W^T) |mv >* ^ mv| 
rn^n v m n 

appi 


(3.60) 


The expression for A g that has been obtained will be studied 
in some detail. In particular, we have 

U 2 ,a 1 ] = [a 2 ,a] = 0 ( 3 . 61 ) 

Thus, the starred basis is constituted by eigenfunctions of A^ 
so that in order to go to higher-order theory another basis is 
unnecessary. Since we can write 

A-Jma)* = 5™|ma>* (3*62) 

we have 

A 1 = Z B mct,mpl ma> <rap l (3 »63) 

m 
a£ 

A 1 - I B m a, J ma> <*l 0.64) 

a? 
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The adjoint statement is 


< mv l A l =^B ma;mp <mv|raci> <m3| 
ap 

0- 6 5). 
P 

which can be written as 

= .< mv l A 1 ^3*66) 

3 

Similarly the matrix representation of the direct statement is 
given by 

A”|rav> =^B majmp |m4> <m3|mv> 
ap 

■ 1 ^ (3 - 67) 

a 

so that 

A“U» = ^ B mamV |ma> (3.68) 

a 

We notice the following relation 


2_j B ma.,nP®nP,nv 
m/n 
op 


v 


lm4> <(nv| 
a -a 


m n 


. V R n lnv> <mai = 

+ nV,nP nPjtna a m “ a n 
m/n 
op 


’X 

msfin i 


B 


ma, nP 


lma> <npjA^ 

a -a 
m n 


,v B 

h L nP^ma 
m/n 
ap 


A^|nP> <ma| 


a -a 
m n 


96 



' X B ma,nP 
m^n 
csp 


|mc£> <nP|A 1 - A x |na> <np| 
( a m “ a n) 


" V "r lma> <np[ ] 

- L-j ma,n 0 ( a m“ a n ) 
m^n 
aP 


(3-69) 


We can then write the first-order equation as 


PtL 


i^l + A.u = rV-R t ma > < np l a lu f t ) 

^t 7 Vl f a -a ) J lj 0 1 

1 m/n V 

ap 


pm^m 


- 1 1 I . B ma,nP B nP,mn?^a^ l mV >* ^ U 0^1 > + 


m/n v^7 

app 


C^o" 1 


- tl I B n»x,nP B pp,ma T^-aJ l raV >* ) 


m/n v^7 
aPp. 


(3.70) 


A slight rearrangement yields 


du, 

+ A 1 U 1 = 



met, np 


|ma> <nP| 

'W 


9^1 


n. 


A> 


-II b 

m/n v/7 
aPp 



(3.71) 


We can then obtain an explicit representation for the first-order 
operator. 
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A iA.,T, a V B mn , v,ft P T i iA-,Xr 1 -iA,X. 

IU-^t-l) = ie 1 1 U 1 ('t 2 )s); -|rg -^j J o dXe 1 [|ma>< n pi ,Aje 1 U q (t 1 )+ 


m^n 

ap 


m/n v^Ay 
iA n t-i 


n C^yC^ p T 1 iA n X „„ — IA-, A. 

, B ma,nP B n|3,mii Ta^^) f Q 6Xa Uv><m 7 |e V T 1> 


-iA, X 


■ tx[e lfllA |-X„ P U-“^j Uo (^ l)+ 

m=/n 

cp 


,m + „m 


m/n v^y 
apjjL 


- 1 ! I V ' 7 I <" v >* tmo-l U Q C=?i ) 

w,_A~ m_/U. in n 0 

(3.72) 


After some rearrangement, we have 






q r IA-,T n -iA n TT 

TC 1 1 ’- a v e 1 | raa><nP| e 1 ma><n3| 

m/n m n 

cp 




+ . /»m fRix 

C m c m / i x 5 v~ 5 'v^ T l 

-Z I A.o,npBnP,^ (3-73) 


m^n v^y 

tpn 


This expression can be manipulated into the form 
IAt t-, v v 1 B. 




m/n 

op 


IAtT, -iA-.T n . 

e 1 1 |ma> <np|e 1 ^(t^) + 


0“> 


V V B B " H. . (e_ 
L L, ma J np'°np,m}J- (a -a ) V , 

m/n v/ 7 m n 1 

cpn 


. t * m *nu 

i(« v -« 7 ) x_ 1 s. 


{«v- 6 ”) y 


JlmV^IjiryiUp^)* 




ma,np 


-i(a -a„)T n 
; ' m n' 0 

(a -a 1 
v ra rr 


lma> <nP|u 0 (T 1 ) 


(3.74) 


We can therefore write. 


B_ 


V*0> - •“ lT \^ a ) + I ^#T I °Sv° 7 P sl (S ”~ 5 ” hl I mv>**<n 7 | V*1 ) + 

m^n m n v 7 
ct£ 


^""1 

4?n 

ap 


-i(a -a )t 

' m r i t 


m n ; On + 




'ma.nP (a -a ) 

3 ' m n ' 

n#n V7 


IC G ?P hvr<n 7 |u 0 (^) 


c m+ c m / WWh, \ 

- I I Ke.r4Prt.rn T^ ^V fe m < '0 1 » v >~<«l V*1 ) 

J3L 77_/„, ' m n /N 1(0 v -o^; 


m/n v^7 
aPjj. 


(3.75) 


We return to the original noninteraction representation operator 
and see 


Vl<V ■ ^ c»; 0?p [ 

wh - ** A * •*• 


> + „n r i ( a m- a n ) T 0 i ( 4 v-® 7 ) 1c i 


-1 


lmv>**<n7|V 0 (i 0 ) + 


ra^n 

ap 

V'y 




4 / pRI ,m N 
G m pm x(o v -o 

°av° 7 ii fe 7 1 -ll 


+ i 


^m^nP^nP^mix 


m/nlyV ' ‘ <V a n )L < 5 v- a 7 ) 

OPil 


mV>^tm7lv 0 (T 0 ) 
(3-76) 


We also have 
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(3-77) 


iAx n iA i' r 1 l*o T c J 

W- # 8 8 W-- 


This gives a compact representation of the uniformized result. 

It is of interest to note the expression for the secularity 
in second order. It is given by 


U 2^ t o) U 0 ^ 0 > _ U 2 f ' T l^ U 0 ( T i B m a,nP B np,nv(" 




(a -a Y 
v m n J 


y|mo><n^ 


V 

^ /„ i(a m - a n )T 0 ,N 

L i / )g ) |nV > <ffla| + 

' m n ' 


-i f a -a )t~ 
/ ' m n' 0 

B ma s n?k • 


•^)|ma)> <nP| + 


1 r* V / '^V a n^ T 0 ,v 

2 Zj La B m / a J npPm / a;n‘'P'\ 'a -a ' J\ a ,-a , 

m/nmVn' m n x m n 

ap a'P' 


m n' 


• |m a> <nPj|m , a7 <n'£'l (3-78) 

The main result (Eq. ( 3 . 77) shows a modified Schwinger-Dyson 
expansion is the propagation operator which is uniformly valid, 
providing double and triple summations that have been considered 
are really convergent. -The main results are explicit formulae 
for the secular parts which can be extracted without ambiguity. 
The main tool employed has been the completeness of the eigen- 
functions of the unperturbed operators. Note that if the starred 
basis still had contained degeneracy relative to A x , then the 
procedure outlined here to remove the degeneracy with could 

be repeated with Ag or higher-order A‘s until the degeneracy 


is removed . 
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CHAPTER V 

TECHNIQUE FOR TREATING 
THE LONG-RANGE KERNEL 
IN THE FLUX EQUATIONS 



The equation for the cosmic ray flux in a random magnetic 
field can be modeled by 



where f(t) is the flux which is considered here to be only a 
function of time. The integral term on the right represents the 
effect of the random field on- the flux; the parameter, e , which 
measures the coupling is considered small. The problem is to- 
determine the long-time (or near eqauilibrium) behavior of the 
flux, which then can be used to determine the transport properties 
of the cosmic ray density. Perturbation or iteration techniques 
produce expansions of f(t) in e which are nonuniform for 
large time, and the standard adiabatic approximation fails because 
of the long range of the kernel which decays only — 1/X for 
large X 

- The type of equation modeled by Eq. (l) appears often in 
studies of kinetic theory, but usually with a kernel, K(X) , 
which is short-ranged due to the typically short-ranged inter- 
actions which are studied. The adiabatic approximation can then 
be made to determine the leading long-time behavior of f(t) 

Since the coupling parameter is small, we assume that f(t) 
changes slowly compared to the time scale over which the kernel, 
which represents the short-ranged interaction, decays. Then a 
good approximation to f(t) is given by 

- - eAf (2) 

where A is a constant given by 



The flux decays exponentially to zero with a rate determined by 
A . Clearly, this approach fails for K(X) = 1/l+X , since A 


101 



becomes, undefined. The problem with A is directly related to 
the decay of K(X) for large X , which is too slow for K(X) 
to be integrable in the sense of Eq. (3). 

The long tail in K(A.) for the cosmic ray flux is due to 
those particles which have very small velocity, compared to a 
typical particle, parallel to the mean magnetic field-. Thus, 
although their velocity may be high, their translation in space 
is quite slow and their interaction with the random magnetic field 
in a given region of space can last arbitrarily long. In previous 
attempts to construct kinetic theories for cosmic rays, the adia- 
batic approximation has been applied, in spite of its failings, 
with artificial cutoffs in the upper limit of Eq. (3) or by 
ignoring the particles in that region of phase space which have 
pitch-angle near 90 ° and which contribute the tail in K(^-) . 

In this chapter, we present several approaches which give 
the long-time behavior of f(t) . One of these approaches, 
which involves the time scale extension technique, is general 
enough to be applied to other kernels which may be more complica- 
ted but which have the feature that for large X , K(X) ~ l/X . 
This approach is ideally suited to the true cosmic ray flux equa- 
tions . 
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1. THE LAPLACE TRANSFORM 


We can obtain an indication of the likely large-time behavior 
of f(t) by studying its Laplace transform near the origin in 
the transform variable, cd . From Eq. (l), we find 


where E-^oi) 


f(a>) 


f (0) 

tt) f •, 

4- ee E 1 (ai) 


is the exponential integral 



(i.i) 


( 1 . 2 ) 


E. (cd) has a branch line extending from the origin along the 

negative real axis to infinity. Hear the origin, it is given by 

00 

%(<!>) = -7 - in co - y (1 arg o>| < tt) (1-3 ) 

n==l 


and for large o> it behaves asymptotically as 






(1.4) 


V is Euler's constant; i.e., y =• .577... . The Laplace trans- 
form, f(ci)) , has a pair of poles very near the origin in the 
negative real half of the oo-plane. We ignore all other poles 
to the left of these, and also the contribution of the branch 
line, to determine the large- time behavior of f (t) . This 

procedure cannot give a proof of the behavior of f(t) unless 
carried to much more detail than we have managed; however, it 
can serve as a guide to the kind of solution we should look for 
with other techniques. 

We look for solutions to 

+ ee ai E 1 (co) = 0 (1*5) 


with 



and with R « 1 . Prom Eq. ( 1 . 3 )* we let E-^cd) = -Jin Tcd where 

r = in 7 • The real and imaginary parts of Eq. (1.5) can be 
manipulated to obtain 

R 2 - e 2 e 2R cos 0 £(in TR) 2 +0 2 J (1.7) 

and 


0 = a + r sin 0 


where 


cos a = 


in TR 


y (in TR) 2 h-0 2 


( 1 . 8 ) 


(1.9) 


For small R and |0|< v , Eqs. (1.7) and (1.8) have a pair 
of solutions which can be written approximately as 


cu = e | in e | e" 

T* 


+ iTTl 1- 


in e 


( 1 . 10 ) 


Thus, for large t , we expect 

f(t) ~ e” e ^ n ^[A cos eirt+B sin eirt] (l.ll) 

where A and B are constants. 

Equation (l.ll) has two surprising features. First, there 
are oscillations on an (et)-time scale and, second, the solution 
is exponentially damped on a el in ej t- time scale. Since the 
damping is fast compared to the oscillations, the oscillations 
can only make a minor modification to the basic exponential 
solution. We can begin to see why the adiabatic solution can 
never give a good approximation to the long-time behavior. In 
that approximation, we effectively impose exponential damping 
on the slow (et)-time scale but the damping coefficient becomes 
infinite . What that solution is trying to tell us is that the damping is 
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finite but on a faster time scale than (et) . All hand-waving 
arguments which attempt to produce a finite damping coefficient on 
the (et) scale actually carry us further from the correct long- 
time behavior since they actually slow down the exponential damping 
even further. Relative to the cosmic ray problem, these attempts 
make the effects of the random field on the particles appear much 
weaker than they actually are . 
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2. THE OUTER EXPANSION 


In this technique , we introduce two tricks which make it 
possible to find the long-time behavior of f(t} . In the first, 
we introduce a new unknown function, g(t) , through 

f (t) = e efc g(t) (2.1) 

Then 

If + €g = -€ J\x g(t-X) (2.2) 

In addition, we introduce 

t = et , g(f) = X(?) (2.3) 

to obtain 

+ x = -e € J +T dX X(e+x-X) (2.4) 

Thus, Eq. (2.4) for x is m°re amenable to the adiabatic approx- 
imation due to the exponential damping which we have forced into 
the kernel. 

For large t , in order to obtain the leading contribution 
to x in an asymptotic expansion in e , we write 

|| + X = - /°V ( 2 *5) 

The adiabatic approximation now gives- 



( 2 . 6 ) 


or 
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(2-7) 


-[1+E (e)l-r 

X(t) - e 


X(0) 


Then 


f(t) = e” €lin €lt f(0) (2.8) 

In this case, the oscillations are missed, but since they are very 
slow in any case, this is not a serious shortcoming of this tech- 
nique. In fact, we will now proceed to refine this technique to 
regain the oscillations; however, it is doubtful whether this re- 
finement could be applied in the case of a more complicated kernel. 
In the next section, we will present a technique which gives both 
the proper damping and the oscillations and which seems more gener- 
ally applicable. 

We go back to Eq. (2.5) and assume the analyticity of x(^) 
Thus, we can Taylor expand xC'T”^) about x(^) • In this case. 


CO 


c>X . 
+ 


1+E 1 (e) X= -) e n a n ” 1 (e) (2.9) 

J J i * OT 


n=l 


where 


a ‘ 1 " 1 (e) = (n-l)le" n 


e T € n_1 

l+e+ 21 • * • + "Cn-ijr 


(2.10) 


We once again consider only the leading terms in e to find 


00 


+ [i +El(e) } . -I 


( 2 . 11 ) 


n-1 


We substitute 


- (l+E n (g)-v)t 

X = e x x(0) 


( 2 . 12 ) 


to determine the correction, v , to the basic ej^n e| decay we 
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found from the adiabatic approximation. .We find 


Thus, 


v = .*n[v-E, (e)3 


v = in | in e | + iir + 0 


0n e 


Neglecting in|in e| compared to E-^(e) , we find 

_ e “(l+|in €|+ i7r)T ^(o) 


and 


f(t) = e - e|in s|t 


A cos eirt+B sin eirt 


] 


in agreement with the Laplace transform evaluation, 


( 2 . 13 ) 

( 2 . 14 ) 

( 2 . 15 ) 

( 2 . 16 ) 
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3- LINEAR TIME-SCALE EXTENSION 


In this section., we present another technique for obtaining 
an asymptotic expansion of f(t) valid for -large t . This 
technique has the advantage of being somewhat more rigorous than 
the approach in Section 2 and also - considerably more adaptable 
to cases where the kernel may be more complicated. 

We introduce an extension of Eq. (l) as follows 


D ^( T 0 J ' r 1 ) 

Dt 



f (Tq— X, T-^— OlX) 

1+71 


(3-1) 


where f is a function of two independent variables, 
and , with the additional condition that f(t,at)=f(t) 

The total derivative, D/Dt , is along a trajectory in the two- 
dimensional (' r o-’ T l^ space specified by 


D _ d , a d 

Dt c)t^ 


(3-2) 


Notice that if we set Tq = t and. = at in Eq. (3*1), we 
regain Eq. (l). This is the sense in which we call Eq. (3*1) an 
extension of Eq. (l). 

The parameter a is meant to be a function of e which is 
small and which is to be chosen to remove the nonuniform behavior 
in time from the perturbation expansion of f . If the kernel 
were short-ranged, we would choose a = e and then the slow time 
scale would be t 1 = et along the "restricted" trajectory. The 
results of uniformizing the perturbation expansion would be 
equivalent to the adiabatic treatment of the original equation 
for f discussed in Section 2. In -the ease of this kernel, we 
find A = co and therefore conclude that we must search for a 
different slow time scale which is, nevertheless, faster than 
et . Thus, we leave a an unspecified function of e which 
must, however, satisfy the following relationships 
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a 


-> 0 


4 



e .j 


(3-3) 


We insert an expansion of f into Eq. (3*1) which is given 
by 


f = ? 0 + I + 


(3.^) 


and equate the coefficients of equal powers of the relevant para- 
meters. In doing this* we must be especially careful since we 
already know that the integral on the right side of Eq. (3«l) is 
not 0(l) but is some other formal order, a , which we have not 
yet determined. Thus., we have 


and 


* f 0 


= 0 


( 3 . 5 ) 


df n e cif. bf 1 

n 0 e 1 , ^ 1 


■-./ 

J 0 


r 0^ f 0 ('v 0 -\,x 1 -aX) g 2 px q f 1 (T 0 -?i J T 1 -aX) 

Q A — — « A 1 > 


1 + x 


a Jr 


1 + X 


( 3 . 6 ) 


We wish to obtain an asymptotic expansion of f to 0(e/a); 
therefore , we neglect the term e(df-^)/dT^ and also the second 
integral term. Then, from Eq. ( 3 . 5 ) 


Sf 0 e Sf l 
a ^F7 + a^T 



1 o^l”^’^ 
1 + x 



f 0 ( T i- aX ) 

1 f 


(3.7) 


Upon integrating over , we find 
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a *l(V T l> “ a + 



1 + X 


4- 


' r o[ a + e f 0 dX 


V-V 0 ^! 
m 


(3-8) 


To remove the nonuniform secular behavior of f^ , we set 




? 0 ^> 
1 + X 


(3*9) 


and, therefore, 

n r ° pco f 0 (r,-cc\) 

f l<V T l> = + a J Q ds J s dX T+X (3.10) 

We determine a from Eq. (3*9) by setting 

f 0 (0,T 1 ) = e Tl f 0 (0,0) (3-11) 


a = -ee a E 1 (a) (3-12) 

The roots of this equation which satisfy the conditions of Eq. 
( 3 .I 3 ) are exactly those which we found in Section 2 in our 
discussion of the Laplace transform of f ; they are given by Eq. 
(1.10). Notice that these roots have negative real parts; there- 
fore, even though f Q grows exponentially with , along the 
restricted trajectory where = ctt , f^t) will decay expon- 
entially as e(J5n e)t . In addition, we find from Eq. (3*10) 


fi(V T i) — ■^i(^ j,r ]_^ 


; i +a 

e 


E 2 (a(l+T 0 ))-E 2 (a) 


(3-13) 


where 



(3-14) 
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For small a } but large Otq , Eq. (3*13) is given approximately 

by 

- T +otr 1 

fl( V Ti) = f^O.Tj) - e - lj (3.15) 

Thus, we have indeed uniformized the expansion of f to order 
1/Un e| . 

Along the restricted trajectory f Uqj ) = f(t) . Thus, 
we have 

Sf 0 

■gx af> 0 ” 0 (3*16) 

and 

f 0 (t) = e afc f o ( 0 ) (3.17) 

where a is given 'by Eq. (1.10). By insisting on the reality of 
f 0 ( t ) , we find 

f(t) - e" € ^ n e ^[A cos eirt+B sin eTrt] (3*18) 

for 


t 


1 

eUn e 


(3.19) 
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4. CONCLUSION 


Of the three methods developed here to obtain the leading 
large time behavior of f(t) , the time-scale extension seems 
the most easily adapted to other more complicated kernels. On 
the other hand, if the slow oscillations were, in some situation, 
considered unimportant, then the outer expansion technique could 
be equally well adapted . 
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CHAPTER VI 


GENERALIZED LARMOR THEOREM 



1 . INTRODUCTION 


✓ 

a . Alfven's Guiding Center 

We split the world velocity of a charged particle in an 
electromagnetic field into a space-like and a time-like orthog- 
onal 4- vectors (with two degrees of freedom each) using field 
strength tensor and its dual. For uniform 6 and B fields, 
the (Lorentz invariant) lengths of each of- these two velocities 
are constants of the motion and generalize the "parallel* 1 and 
"perpendicular" energies in a constant magnetic field. The 
space-like velocity always contains the periodic (gyratory) part 
of the motion of the particie, while the time-like component is 
the world velocity of the particle's generalized guiding center. 

The equations of motion for the two projected velocities 
are solved in a " standard " configuration and combined to give 
the general solution for the trajectory of a particle in an 
arbitrary uniform field configuration, with the help of a simply 
constructed Lorentz transformation that represents a general 
e x B drift. 

Alfven's definition of a guiding center description of the 
behavior of a charged particle (relativistic, or not) in a uni- 
form magnetic field follows naturally from basic properties of 
the particle's trajectory. Specifically, the magnitude of the 
velocity is conserved in time and the motion is easily visualized 
in terms of gyrations about a point which moves with constant 
velocity along the magnetic field. The cycle average position 
of the particle is well-defined and called the guiding center. 

The distance from the guiding center to the particle position v 
is called the gyration radius. In situations where, in addition 
to the magnetic field, there are weak electric fields and/or 
small gradients and/or time variations in the magnetic field, 
the Alfven^ guiding center plus gyration vector approximation 
to the particle motion is valid. In this approximation, the 
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guiding center retains its usefulness; the first-order correc- 
tions to its motion are computed using standard perturbation 
theory in terms of the small field variations. Even to first 
order, this technique has been very useful in many applications 
(see, for example, Northrop and Teller 2 for a study of the 
trapped particles in the Earth’s geomagnetic field). More 
recently, higher-order terms than the first have been calculated 
by Kruskal^ for a general class of dynamical systems which 
include the charged particle in a field. 

b. New Results 

In this chapter, we re-examine the study of the guiding 
center in uniform and stationary electric and magnetic fields of 
arbitrary magnitudes and directions . We obtain a relativistic 
generalization of the guiding center in space-time which is 
identical, in the nonrelativistic limit in 3-spa ce * to the Alfven 
guiding center. 

We employ an operator formalism using projection operators 
which are constructed out of the field strength tensor and its 
dual. The pair of projection operators which we use form a 
complete orthogonal set which allows us to project the world 
velocity of the particle into two parts which lie in orthogonal, 
2-dimensional subspaces. One of these parts is the world velo- 
city of the guiding center and is time-like. Thus, in any non- 
singular field configuration, we avoid the problem of averaging 
over the particle trajectory to obtain the guiding center posi- 
tion by merely projecting away the oscillatory part of the par- 
ticle’s velocity; the part that is left represents the (time-like) 
4- velocity of the guiding center. 

The operator formalism constructed in this report is applic- 
able to any uniform field configuration except the singular case, 
S*B - 0 and B - 6 = 0 , where the operators become undefined. 
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Notation 


c . 

Gaussian units are used throughout this paper. The trans- 
posed matrix is denoted by a tilde and a prime is used for a 
Lorentz-transformed quantity. The determinant of the matrix A 
is written det A . The dual of a tensor, defined by Eq. (1.3) 
is denoted by a star. The world velocity, u , is written as 

u = ( 7 v,ic 7 ) , 7 = [l - v 2 /c 2 ]“ 1//2 (l.l) 


and the proper time t is related to laboratory time by 

7dt = dT (1.2) 

We define the dual f* of a skew tensor f as 


f* == ± € f . f** _ f 

liv 2 |i.vaT] arj ' x 1 


The dual tensor is sometimes defined as 


F*;=(i/2i) Vxp f Xp 


Then 


(f*')*'= -f 

A 

A 3-vector is denoted by A and its magnitude by 

A = (A- A) 1 / 2 

\ 

The signum function, sgn (x) , is defined by 


sgn 


r +1 x > 0 

< X > =1-1 X < 0 


(1-3) 

(1.4) 

(1.5) 

( 1 . 6 ) 

(1.7) 


116 



2. ELEMENTARY PROJECTION OPERATOR FORMALISM 

We use the motion of a nonrelativistic charged particle 

in a uniform magnetic field to motivate the general projection 

operator formalism presented in the next section. Standard 

analysis is reproduced with (a special case of) our technique 

which projects the behavior of the guiding center and of the 

gyration vector from the equations of motion and thus facilitates 

carrying out their solution. To this end., we write the equation 

of motion in terms of (a 3 x 3 submatrix of) the electromagnetic 

field tensor (see Eq. (3-2)). The properties of the particle 

motion are then characterized through the algebraic properties 

i 

of the frequency matrix. The velocity vector, v , satisfies 


d^ v = v , 


m Q c r B 3 ° 

\ B 2 ~ B 1 


( 2 . 1 ) 


The useful properties of the field strength tensor are: 
(A) The matrix O is antisymmetric . Thus, 

0 = -0 


( 2 . 2 ) 


and therefore Tr 0 = det 0=0 . We also have the represen- 

tation 


n i.1 - e i.1k“B P k ’ ^k " V 6 » “B - 


(2.3) 


Furthermore, exp (Ot) is an orthogonal matrix which represents 

~ a 

a rotation about the direction p with rotation angle (Bt) 

exp (fit) = R(p,Bt) = R((3,-Bt) (2.4) 

(B) The matrix 0 . 0/4^ is a projection . Thus, 

Hs - = I - PP , o£ = 1 Tr 0 • O (2.5) 
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is a projection operator (normal to the magnetic field) since 
it is idempotent and symmetric 

N • N = N , N « N (2.6) 

It follows that N is positive definite 

N > 0 (2.7) 

and that 

P = I - N = 6P (2.8) 

is a projection that satisfies the commutativity and orthogonality 
rules 

N • P = P • N = 0 (2.9) 


(c) The matrix 0 is "parallel" to N and "orthogonal! 1 to 
P . Thus, 

N ♦ 0 = ^ • N = ^ , P • 0 = 0 • P = 0 (2.10) 


We can now show that the basic properties, of the particle 
motion reflect our algebraic rules. For this purpose, we use 
the two projection operators introduced, N and P , to decompose 
the velocity vector into two mutually orthogonal parts 

v_j_ = N • v , v |( = v - v ± = P ♦ v (2.11) 

each satisfying Eq. (,2.l) because of (2.10). 


We obtain the standard particle properties in parallel 
with the matrix properties just- discussed. 

O . 

(A) Conservation of kinetic energy. The quantities v , 

Q T-— — 

V 1 * \\ are conserve(3 because of (2.2), e.g.. 


1 d , 2 \ 

2 3? ) = v 


dv 

dt 


a 

= v 


a 


a 

v 


- 0 


(2.12) 


Furthermore, the noninert ial effect of the constant magnetic 
field is eliminated by transforming to a reference frame rotating 
with angular frequency oig (Larmor theorem). This immediate 
consequence of (2.4) is generalized in the next section to include 


118 



relativistic motion as well as electric fields. 

(B) Periodicity -of Vj_ with angular frequency tOg . Since 
satisfies (2.1), we have, using the projection property {2-5), 


d 2 a ^ d a ^ _ a 2 a pa 

^5 v ± ■ S- « v x = § • 5 • V X = -®B * • v x = -®B V X 


(2-13) 


The vector v ± is harmonic because of (2.7). 

(C) Inertial motion of the guiding center. Using (2.10) 


d u „ i 

at v ii = P • e • V = o 


(2.14) 


J A 

Thus, the integral of v^Ax^ represents inertial motion. 


Integration of 


a a -i ^ AA-i d A 

V = v+ - 2 n.n.v = v + - 5 fi. 3F v 

“b ®b 

motivates the definition of the guiding center by 

A A 
A A % X V 


(2.15) 


.A A I _ A A 

Ax (] = X H 5 Cl • V = X - 

“b ‘ 


“b 


2— S A *GC 


(2.16) 


The magnitude of the gyration vector, x s Ax ± , is the usual 


radius of gyration r and is given by 

s 


r g - l x ±l = v i< 




(2.17) 


Combining the integration of (2.13) and (2.14), we obtain 
useful representations for the particle trajectory 

v(t) = (p + N cos oigt + ~ £ sin oigt).v 0 (2.18) 


Ax(t) = v |I0 t + 


sin cDgt 


cos 


m 


I - 


v- : 


B 


<4 


£2 )*v 


±0 


(2.19) 


The projections P and N can be eliminated in favor of Cl by 
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using (2.8) and (2.5)* In the next sect ion we show that the 
inclusion of relativity and electric fields is obtained by a 
straightforward analogy. 



3* GENERAL FORMALISM 


The Lorentz equations of motion of a relativistic charged 
particle in an electromagnetic field are 


QT 


U = 03 • U 


(3.1) 


where 03 is the antisymmetric field strength tensor 


/ 0 =3 - B 2 - 1S l\ 


0 ) - 


m 0 o 


-B 3 0 


B 2 - Bl 0 -i® 3 

i® 2 18, 0 


= -CO 


1 


(3-2) 


From Eq. (1.3)* to* is given by the substitutions <-*■ ~iS^ 

From to and o>* we can construct two invariants and two useful 
tensor relations 


i=* a. ^ O 

det 03 = det co* = -( — ) (g.B)^ 


m 0 c 


2 1 

03 = i 03 03 „ 

2 [xv |xv 


= ~ 03* ,03* , = - g 2 ) 

2 |XV [XV Vjn n c/ V ' 


03 * 03 + 03* *03* = -03 I 


(3-3) 

(3.4) 

(3.5) 


A ^ ^ 

03 *. 03 = 03*03* = i(g.B)(j~-J I = [ sgn(g*B)3 (det 03)1 (3*6) 

Be low , we obtain two orthogonal projection operators analo- 
gous to those of the previous section 


EL + n l i = 1 - IL-EL =Hi -IL - 0 


(3-7) 


which separate u into two orthogonal parts 


U X 3 IL * u - u n = 


(1 - u - Ujl 


(3-8) 
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0 


(3-9) 


u. 


tt, = u, 


u„ = 


in such a way that u ( j represents the generalization of the 
standard guiding center velocity. The instructive special case 

-X -A 

£ * B = 0 is discussed in Appendix F. 

Consider the operator 



- 0>CD+Xl 

(XT 

1 + 2X 


where X is a scalar determined from the condition 

ifc®> • lit®*) = lie®*) • 11(5) - o 


( 3 - 10 ) 


(3.11) 


Notice that we have normalized so that 

]](“) + ]fe*) = 1 (3-12) 

Prom Eqs. (3.10) and (3. 11) we find that X satisfies a quad- 
ratic . We choose the root 


X = 






1+4! 


J. A 

g»B 

■b 2 -s 2 ' 


(3.13) 


a i 

to satisfy the condition that X ~> 0 as g*B -* 0 since in this 
case (3*6) yields the projection property, i.e., (3*7). Using 
Eq. (3*13)., we obtain 




(3.14) 
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The useful matrix properties of E(9) are given in Appendix 




G. We note that if fi-B ^ 0 


Lim = Lim J](cp*) 


B 2 -e 2 -o + 


B 2 -£ 2 ~<r 


(3-15) 


Lim JJ(ap) = Lim tt«*> 


b 2 -£ 2 -*o" 


B 2 -e 2 ~>o + 


(3.16) 


Therefore, we define the projection operators, IL and Ei * 


IL- 


11(9 )B 2 > £ 2 

E(9 *)b 2 < e 2 


IL « 


11(9* )B 2 > £ 2 


E(9)b : 


2 < g 2 


(3.17) 


Then El and Ei are continuous at B -£=0 and they have 

the desired properties defined by Eqs. ( 3 . 10). Prom (3*l); we 

see that u ± and u (! obey the same equation of motion as 

u(x) does. Therefore, if Uj_(0)=0 , then u_ l (t)=0 for all x . 

Thus, the guiding center represents a particle for which u JL (0)=0. 

Since u. is space-like (see (4.10)), it cannot be directly 

2 2 ' 

related to the motion of a particle. Furthermore, u (| and u. ± 
are Lorentz and gauge-invariant constants of the motion. 

By differentiating Eq. (3 • 1 ) 


d 2 u, 


d 2 u, 


dx 


CD . o> • u x , 


dx 


= CD ♦ CD • U,, 

2 II 


(3.18) 


But, from Eq. ( 3 . 10) 

co • cd = cd 2 [XI - (l-!-2X)J|(cp)] (3.19) 
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Therefore 


d2u ± 2 1 

57“ = v 


d 2 U|, 2 

— - % U I! 


a l 


dT 


«5 - 


CD 2 (l+A.) 

B 2 

> e 

cu 2 X 


> B 

o> 2 X 

B 2 

> e 

o) 2 (l+X) 

. £ 2 

> B 


( 3 . 20 ) 


( 3 - 21 ) 


The solutions of Eqs. (3*20) and (3«2l) are 


Uj_(t) = 

u„(t:) = 


sin CD.T + I cos 03* T 


03 

03 ,. 


}IL- u (°) 


n 

sinh 4- I cosh co^ tt J *H^| * u(0) 


Therefore 


CD 


± 


( 3 - 22 ) 

(3-23) 


e CUT = sin cd^t + I cos + s i n ^ CJ u T + ^ cos ^ tu n T )W\\ 


% 

(3-24) 


Integrating the resulting expression for the velocity, we obtain 


r /sin cb. t\ 

- -(o) = K-^rv 


CD* 


v'deFTo 


(cos cb ± t - 1) *Uj_(0) + 


'sinh <0,,^ + 


(B 


II 


Vdet co 


(cosh <b (( t - 1) 


■u„(0) (3-25) 


Below, we use this expression to study the motion of the general- 
ized guiding center.. 
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4. STANDARD AND SINGULAR FIELD CONFIGURATIONS 


The possible field configurations are given in Table 1 . 
Below, we show that only two configurations, the "standard" and 
the "singular," need separate treatment, all others being related 
to the standard configuration either as special cases or by an 
explicitly given Lorentz transformation. 



Table 1. Possible Field Configurations 


(a) The operators constructed in Section 3 are singular ; 
thus, we call this the "singular" configuration. We give the 
solution of the particle equation of motion for this configura- 
tion separately. 

(b) This is the familiar crossed-field case where the 
magnetic field dominates.- The standard technique which is used 
to obtain solutions of the particle equation of motion in this 
case is to consider the motion in a Lorentz frame in which 

= o . The motion in the original frame is obtained by merely 
§dding the constant transformation velocity to the solution: 
u = (c/B 2 )(e xJ3) . Thus, by studying the motion of the particle 
in a uniform B field, we actually study the entire class of 
field configurations. 

A 

(c ) Using similar reasoning, we choose the pure £-field as 
the prototype for this class. 
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(d,e,f) We vfill see that when £ and B are parallel, Jin 
ill reduce to particularly simple diagonal forms which are inde- 
pendent of the magnitudes of the field strengths. Furthermore, 
the parallel field configuration can be reached from any nonper- 
pendicular configuration through a Lorentz transformation given 
in Appendix B . Therefore, if we adopt the parallel field case as 
the prototype of each of classes d, e, and f , we may study 
all three of these in one example. Furthermore, since B 2 - £ 2 
can take on any value in this example, we can also study the 
special cases, B2 = 0 or £2 - o , but these are the pure field 
cases which we have adopted as the prototypes of b and c 
Thus, the prototypes of all classes except case a can be studied 
through one example in which the JJ* s take on very simple forms. 

We now give explicit solutions of the equations of motion 
in the standard and in the singular configurations. 


a . Standard Configuration 

A 

In the standard field configuration (S.C.) we let the £- 
and B-fields be parallel and in the z direction. Then JJx and 
!!,[ are given by 




(S.C.) (4.1) 


Thus, 

u i. = ('YV 1 ,7V 2 ,0,0) , Uj, = (0,0,7v 3 ;i7c) , (S.C.) (4.2) 

i. a 

In the standard configuration, u (I = 7V 1( represents the motion 
of the particle along the magnetic field. Henceforth, we will 
call u (| the world velocity of the generalized guiding center. 

From Eq. (4.2) we see that 

u l> 0 (4.3) 

2 2 2 

but we always have -c = u^ + u I( and therefore 

^ “° 2 ( 4 - 4 ) 

Thus, u [( is time- like in the standard configuration and u ± is 
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space-like. By the discussion of the previous subsection, these 
properties hold also in any configuration. Prom Bqs. (4.2), if 
x(t) is the 4-space position of the particle, then the position 
of the guiding center is given by 

x„ = (0,0,x 3 (t),x^(t:)) (4.5) 

and the position of the guiding center "complement” is given by 


x ± = (x 1 (t),x 2 (t),0,0) (4.6) 

With the fields in the x^ direction, the vector components of 
Eq. (3*25) are given by 


/■sin oi_t\ /cos ov/t-l - ' 

Ax l = x l( T ) “ x l(°) = u l^°\ (3^ J'~ u 2^°\ ", 


(4.7) 


/sin ov-tn /cos 

Ax 2 = x 2 (t) - x 2 (0) = u 2 (0)(— ^ J + u x (0)Q — J (4.8) 


xsxnn /cosh cd«t~±\ 

Ax 3 = x 3 (t) - x 3 (0) = u 3 (0)(^ ■£— j + C7(0)(^— J(4»9) 

o o 

/sinh cd r\ /u 0 (0)\/cosh cd_t-1\ 

* ■ *K— 57 + (-V- X— ST 6 ) (4 • 10) 


where 



(4.11) 


and where we have set- t(0) = 0 . This result is valid for any 

S and B . In addition, without a loss of generality, we may 
set u 3 (0) = 0 , since a Lorentz transformation 'parallel to the 
fields does not affect the fields at all. 


In order to write Ax n 


A x r 


j <— l-A-g j uixv 

it is convenient to introduce the new' variables 


and Ax„ as functions of t 


x (+) (r) 5 Ax 3 + iAx^ 


x<’>(t) = Ax 3 - iAx^ 


(4.12) 


127 



(4.13) 


where Ax^ - icfe . Then 

u (+) = + iu^ , u (_) = 

Now* Eqs. (4.12) are equivalent to 



Equations (4.14) may be solved for cd^t 





(4.14) 


We obtain 


>] 

4 


(4.15) 


Upon eliminating cc^t } we find an implicit relationship between 

Ax~ and t which can be inverted to 
3 


Ax 3 (t) 



(4.16) 


Equation (4.16) is a general statement of the 3“P° s ition of the 
guiding center in the standard configuration as a function of 
time . 

With Eq. (4.16) and x^”^ = Ax^ + ct , we may obtain an 
explicit expression for r(t) 





(4.17) 


This last expression may be substituted directly into Eqs. (4.10) 
and (4.11) to obtain Ax-^(t) and AXg(t) 
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For (cD g t/7(0)) « 1 , Eq. (4.1?) is given by 


OigT - 





'in 




1 + 


“ 

7 ( 0 ) _ 


CD t 

B 

tCOT 


(4.18) 


With this expression, we find 


Ax-^ ~ v 1 (0) 


sin (m B /7(0))t 


- v. 


/cos K/y(0))t - In 

: (0) ( <VWJ ) 


(4.19) 


Ax 2 - v 2 


/sin (civ/yCOjJtx /oos (<rt/7(0))fc - 1 \ 

(°)( 5577{^) ) + V°\ a^TvIOl ) 


and from Eq. (4.16) 



(4.20) 


Note that 



(4.21) 


Thus, the condition that (cD g t/7(0)) « 1 is equivalent to the 
statement that the motion of the particle along the fields is 
subrelativistic . Either t is small enough or £ weak enough 
such that the particle is not appreciably accelerated along the 
fields. The motion perpendicular to the fields need not be 
subrelativistic, however. 

a 

Equations (4.12) and (4.13) are exact when S and aig 
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are equal fco zero. The nonrelativistic gyration frequency is 
corrected to include the relativistic mass; i.e., the relativistic 
gyration frequency is 


_ eB _ eB / 2 |_ 22 

7 ” THIqC “me ' 

If S ^ 0 3 the particle is initially accelerated uniformly 
along the electric field and the motion perpendicular to the 
fields is unaffected. 


When (a>gt/v(0) ) » 1 , Eqs . (4.1) and (4.17) are given by 

Ax 3 r ct (4.23) 

and 


V - 


7(0) 




Thus., in 3“ s P ace * the guiding center moves along the electric 
field with v^ -sr c and the guiding center complement still 
gyrates about the guiding center but at a much slower rate than 
in the subrelativistic regime. 


b. Explicit Tims Dependence for the Singular Configuration 


Let the field components be given by 

e = (B,0,0) and B = (0,B,0) (4.25) 

The components of Eq. (3*1) are given by 

h U 1 - - icQ B u 4 * h U 2 = 0 ^- 26 > 

it u 3 - “b u 1 'ft u 4 = ia> B u l ( 4 - 27) 

For initial conditions, we choose without loss of generality 

t ( t=0) = 0 (4.28) 
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(4.29) 


v 1 ( 0 ) = v 2 ( 0 ) = v 3 ( 0 ) = 0 ' 
With these initial conditions, we find 


£ 

O 

II 

1 — I 

3 

. u 2 = 0 

(4.30) 

u 3 = ( 1 / 2 )c(cq b t ) 2 

, u^ = ic + (l/ 2 )c(oi B T ) 2 

(4.31) 

From Eq. (4.31) for u^ = 

7 = 

= 17 c , we find 
1 + (1/2) (o^xf 

(4-32) 


Thus,, the equation dt/dr = 7 can be integrated directly to 
find the cubic equation 


t 


1 + 




Equation (4.33) can inverted to find 


(4-33) 


(t) = (S^t) 1 / 3 ! (l+ 


8 


(3 V)‘ 


\!/2 j 

1 

0 + 1 J 

T 


1 + 


8 


(3^t)‘ 


1/2 ^l/3\ 


J 

(4.34) 


Equations (4.30) and (4.31) can be rewritten for the three 
velocities. We obtain 


v. 


V 


1 + 


!<‘v) s 



la = 

c 


V 


V. 


(4.35) 


A. 

Equations (4.34) and (4.35) give the solution for v(t) . 

For oigt , Eqs. (4-34) and (4.35) reduce to the nonrela- 
tivistic result 


T~t , ' /"l (m B t)2 (4 - 36) 

•The particle is initially accelerated uniformly along the 
electric field and slowly gains a velocity component in the 
z-direction orthogonal to both the fields. After the particle 
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is accelerated for a long enough time, we have » 1 

Then 


OigTT (ecDgt) 1 / 3 


and 



c 


2 

(6o^tp73 



c 


( 4 . 37 ) 


( 4 . 38 ) 
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5 . CONCLUSION 


The operator technique developed here has been shown to 
be very useful in describing the relativistic guiding center and 
also the general motion of a charged particle in an arbitrary 
uniform field configuration. The relative simplicity of the 
technique leads us to believe that it may also prove- useful when 
the fields are not uniform and stationary. We shall discuss the 
application of our projection operators to the motion of a particle 
in fields which have small variations about the uniform configura- 
tions elsewhere , Since the positions of the guiding center and 
its complement are given in terms of the particle's position by 
simple contact transformations in space-time., it is not at all 
clear that the usual perturbation techniques will be necessary 5 
exact canonical equations of motion for the guiding center and its 
complement may be obtainable , 

Our main results are summarized in Table 2 where the pure 
magnetic field is compared with the general case. 
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Table 2. Comparison of the Pure Magnetic Field with the General Case 
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APPENDIX F 


Note that 


s 1 -- 


0) ■ CD 
2 5 


= 


* 

0) 


CO 


(FI) 


cd or 

form a complete but nonorthogonal set since, from Eq. (3.9) j 

Sr 2 -(T 1 -ST 2 =^i (F2) 


CD 


^2 

Thus w and (# play the role of the projection operators if 
, Furthermore, from Appendix G, 


e i b = o 


c-2 * * 

(A) • CD = CD 


y 1 • 

Introducing 
we find 


CD = cu 


7f 2 • 
Sr 1 • 


CD = 0 , 

(e • 

B = 

0 ) 

(F3) 

* 

CD = 0 , 

(£ • 

_5 

B = 

01 

(F4) 


u 1 => 1 


U 


u 2 = ^ 2 


u 


.du 


dr 


du 

dr 


du 

dr 


= 0 


(F5) 


(F6) 


The vector u^ represents the constant velocity of the rela- 


tivistic guiding center when £ • B - 0 


In fact, the expressions 


u 2 = u - u 1 = u -^' 1 . u = u+ ^ ^).u = u + cd . (FT) 

integrates to’ 

X* = x + -f C (F8) 

CD^ 

.1 

whose 3-part with £ = 0 coincides with the position of the 
guiding center (2,14), but £ . B = 0 
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APPENDIX G. USEFUL MATRIX PROPERTIES 


Using (3.11) and (3.12), we find 

- ft 1 • P 1 -Jr 1 - x , P s . Sr 2 = Jr 2 - ! 

or m 

Also, notice that 

^2 * * ^2 * , /det cd 

(j) • CD — CO * W ss CD + C / rr — CD 

v CD 4 

2 <V"2 /det cd * 

^ * Cl) — CD • --Of tt— O) 

V aT 

and 

& 1 . CD= (0 -t? 1 = »+ 0/^5 <0* 

V QJ 

vl * . * ^1 _ /det cd 

ZT • CO = CD . id = - a / TT— CD 

v m* 


In Eqs. (Gl) and (G2) we have introduced the notation 

a = sgn ( S . b [B 2 - e 2 ]) 

With these relationships, we may go in a straightforward 
to obtain 


If 

If 

IP 


m = 0) . If- 


/det cd „* , 

“ a J Zj- ® + Xa> 

V CD 


CD* = CD* • ]J 2 


• CD = CD • I} 1 = 


1 

+ 2\ 

(1+X)CD* 

, _ /’det cd _ 
+ a/ jp <o 

V CD 

1 

4- 2A. 

(1+X)a> + v /jM ir ®* 


CD 


(Gl) 


(G2) 


(G3) 


(G4) 

manner 

(05) 


(06) 


(07) 
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1 + 2 \ 



and 




/det co * 

- V + xo) 

0 ) 


1 +■ 2X 


(g8) 


We shall use Eqs. (G5) through (g 8) in our construction of the 
equations of motion of the projected world velocities. 
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APPENDIX H 


Given an arbitrary field configuration (g,B 3 0) where 6 
is the angle between £ and B , we can find the equivalent 
standard configuration (£ ' 9 b' 3 0' = 0) . The magnitudes of the 
standard fields are 


/ 


B' 


i(B^-£^)|"l + >/— -(.A + 4(£B) cos 2 q ] 

2 L B 2 - £ 2 


(HI) 


= \{ e 2 - b 2 


)[i + 




/ 2 2. 2 
(b - e 2 ) 


+ 4(£B) 2 cos 2 6 


- B e 


(H2) 


The velocity which accomplishes the transformation from the arbi- 
trary to the standard configuration is given by 

V = a( £ X B)c (H3) 

where 

a = (B 2 +£ 2 ) - J (B 2 -£ 2 ) 2 + 4(£B) 2 cos 2 0 (H4) 

2(£B) 2 sin 2 8 

Thus., for arbitrary (gjB } Q) it is possible using the proper 
Lorentz transformation to obtain ' 3 0) with £ ' and B in 

the z '-direction. Then Hj. an d Hu take on simple forms which 

are independent of the magnitudes £ ' and B' 
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